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CAN THERE BE A SHIELD FOR GRAVITATION? * 


BY 


W. F. G. SWANN! 


I believe the high spot of interest of this organization has been the 
search for a ‘‘Shield for Gravitation.”’ The origin of an expectation for 
a shield is obvious. We can shield our eyes from the sunlight by dark 
glasses. We can shield ourselves from the wind by a suitably placed 
screen. We can shield electrical apparatus from the effects of external 
electrical disturbances. By placing our portable radio inside a metal 
container, we can shield it so effectively from such external disturbances 
that it will not work at all. It will not even respond to the external 
influences which we wish it to respond to. 

Our office is noisy. Every sound we make comes back to us by 
reflection from the walls, so we cover the walls with absorbing curtains 
or felt, and echoes are eliminated. 

In an enlarged sense, we can shield our noses from unpleasant odors, 
from poisonous gases, or from the pollen which carries hay fever, by 
covering our nostrils with a suitable filtering appliance. We can shield 
our heads from the sun by wearing a hat. 

And so, in the light of these experiences, we can well wonder whether 
it would be possible to reduce, or completely obliterate, the force of 
gravity by some kind of a shield. A falling body increases its velocity 
continually under the action of gravitation. Why should it not be 
possible to design, from some suitable material, a container whose walls 


* This communication represents an amplification of an essay of the same title, submitted 
to the Gravity Research Foundation. The essay won one of the prizes offered by the Founda- 
tion and was presented verbally in the form here given at the Foundation’s meeting held on 
August 27, 1960, at Gravity Village, New Boston, N. H. 

1 Director Emeritus, Bartol Research Foundation of The Franklin Institute, Swarthmore, 
Pa.; and a Senior Staff Advisor for The Franklin Institute Laboratories for Research and 
Development. 


(Note—The Franklin Institute is not responsible for the statements and opinions advanced by contributors in 
the JouRNAL.) 
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would ‘‘absorb”’ the gravitational pull, so that the body inside would 
no longer fall ? 

Of course, during the last few years, we have become familiar with 
effects simulating rather closely those which would be caused by the 
reduction or complete elimination of the force of gravity. An animal 
traveling in one of the ‘‘sputniks’’ which encircle the earth experiences 
no sense of weight. If you put him in the center of his cabin without 
support, he would stay there, encircling the earth like a little satellite, 
quite independently of that other satellite which is his container. 
The phenomenon of weightlessness here presented to us is, however, 
not one involving a shield for gravitation, but one in which the animal 
in the sputnik is exposed to two oppositely directed influences, the force 
of gravity, and the effects resulting from the continual change in direc- 
tion which he experiences as he travels around the earth. Indeed, a 
person who jumps off a stool experiences a sense of weightlessness as he 
falls, but it would not be correct to say that the force of gravity has 
been eliminated in the fall. We must rather hold that the force of 
gravity is present in its full intensity, but the acceleration he experiences 
during his fall produces upon all that responds to ‘‘force’’ within him 
an effect just equal and opposite to the effect produced by the force of 
gravitation, so that his condition is the same as if there were no force 
of gravitation and no acceleration. 

If, indeed, we could find a means of shielding off gravitation, or even 
of reducing it, the effects on our everyday activities would be consider- 
able. If we could reduce gravitation to the value it has on the moon, 
for example, a man could jump six times as high as he can at present. 
The limit of his horizontal jump, instead of being about 18 feet as it is 
at present, would be about 100 feet. 

Now when, with the foregoing situation in mind, we approach a 
conventional man of science and ask him to produce a shield for gravita- 
tion, we find that he looks sad and shakes his head. If we ask him 
whether these unhappy and negative gestures imply that he ‘‘knows”’ 
that there can be no shield for gravitation, he will be compelled to 
admit that nobody can make such an assertion. If we ask him: ‘“‘Why, 
then, do you look so miserable ?’’ he will probably be driven to reply 
thus: 

“My friend, in the development which science has experienced during 
the last few hundred years, I and my ancestors in science have learned 
certain things about the way in which the forces of nature act—the 
way in which electricity acts—the way in which gravitation acts, and 
the result of all this knowledge is such as to make me feel that a shield 
for gravitation would represent a very unnatural possibility. It is as 
though you asked me to produce, by cross fertilization, an animal which 
ran on wheels, asserting that such an achievement should be possible 
because automobiles run on wheels. I should have to tell you that the 
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whole structure of animals is such that I cannot see how a wheel can 
ever become a natural part of an animal’s equipment. My knowledge 
of the necessity for blood circulation, and the difficulty of imagining 
any such process operative in the case of a wheel, puts that device into 
a category entirely different from those for which provision has been 
made in biology.” 

And if our man of science, inebriated with the exuberance of his own 
verbosity in giving us this answer, now looks a little less sad, and a little 
more companionable, it may be the time to ask him just what there is 
about the ordinary activities of nature which can look with equanimity 
upon the existence of a shield for certain forces, such as electrical forces, 
and yet compel us to look askance at a shield for gravitation. Having 
asked this question, let us make an attempt to answer it. 

Now let me say at the outset that there has been a theory of gravita- 
tion which could, in principle, lend itself to the possibility of a shield for, 
or perhaps we had better say, an absorber of gravitational force. This 
theory of gravitation, put forward by Le Sage many years ago, postu- 
lated that all space was filled with small particles traveling in all direc- 
tions with high speeds. Let us now consider two bodies, A and B, 


Fig. 1. 


Fic. 1. 


If A were alone, it would be bombarded by particles from all sides 
and would show no tendency to move. This would be true whether the 
particles falling upon it were totally or only partially absorbed. Now 
when B is present, some of the particles which would have come to A 
from the right are absorbed by B, so that A is bombarded more from the 
left than from the right and moves to the right. By analogous argu- 
ment, A partially shields B against bombardment from the left, so that 
B tends to move to the left. The net result is that A and B move 
towards each other as though they attracted each other. If the bodies 
A and B were very thick in the left-right dimension, or if they were 
composed of some material very opaque to the particles, so that indeed 
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no particles reached the space between them, a body placed between 
them would show very little or no tendency to move towards the left 
or right as a result of the presence of external bodies to the left or right 
of both A and B, and we should have a “‘shield”’ for gravitation. 

A hundred years ago, we might have accepted such a speculation as 
is here involved as something not unreasonable. For we recall the fact 
that light radiation from the sky is stopped completely by a very thin 
sheet of metal foil, whereas it can penetrate with ease the whole atmos- 
phere which, as far as mass is concerned, represents the equivalent of 
14 pounds per square inch, and we might well wonder whether there 
might be some substance so opaque to the flying particles that a rela- 
tively thin sheet of it would afford a shield. Alas for such speculations, 
however, we have learned so much about the structure of matter during 
the last century that we see no place for the existence of such a sub- 


stance. Neither do we see any place for the flying particles. Our 
studies of cosmic rays have indeed introduced us to real particles of 
enormous energy per unit mass flying through space in all directions. 
However, our critical faculties have become so sharpened in the studies 
of these rays that, today, much more than formerly, we must look with 
considerable scepticism on the existence of such particles as would have 
the properties necessary to account for gravitation in the foregoing sense. 

Let us, however, turn to the crucial differences between electric 
forces, on the one hand, and gravitational forces on the other hand, 
which make it so easy to understand shielding in the former realm and 
so difficult to understand it in the latter. 

If A, Fig. 2, is a positively charged body, S a metallic shield, and O 
a point within it, we know that apparatus at O will be “shielded” from 
the influence of A. However, we do not visualize the shielding pro- 
cedure as one in which the electrical forces emanating from A are some- 
how or other absorbed by S. We recognize the known fact of the exist- 
ence of two opposite kinds of electricity, positive and negative. The 
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shield contains both kinds in equal amount, so that it is neutral as a 
whole. The positive electricity is fixed in the atoms of the “matter” 
of the shield, but in an electrical conductor the negative electrons are 
free to move, and they do so move as far as they can, toward A which 
attracts them, thus leaving a positive charge on the part of the shield 
remote from A. Now if we should put a charge of positive electricity 
at O, the force due to A would push it to the right. However, the push 
due to the charge labelled Q, and the pull due to the charge labelled Q_ 
combine to urge it to the left, It turns out that the push to the right 
resulting from A just balances the pushes and pulls due to Q, and Q_, 
respectively, and so the charge at O is “‘shielded’’ from the external 
influence originating at A: We see that in this case the ultimate 
phenomenon of shielding depends upon the existence of two opposite 
kinds of electricity. If there were two kinds of gravitational matter, 
one kind K, which repelled its kind, and another kind K:2, which re- 
pelled its kind, and if the kind K, attracted the kind K2, then we might 
hope to realize shields for gravitation in the same manner as we realize 
shields for electrical forces. 

Consider another rather more subtle case, where a beam of light 
waves or radio waves impinges upon a plate P, Fig. 3, with the result 
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that nothing is observed to the right of P. The accepted way of regard- 
ing this phenomenon is not to suppose that the waves are somehow or 
other absorbed by P. We, indeed, say that nothing happens to the 
waves. They go on right through P, off to infinity. If we ask: ‘‘Why 
then do we not observe them?,” the answer is that they set up electric 
currents in the plate P, which currents send out waves on their own to 
the right of P, but waves which are “out of step’’ with the original set 
in the sense that the ‘‘crests’’ of one set coincide with the “troughs” of 
the other. Thus, while there is ‘‘nothing”’ to the right of P, the nothing 
is, as it were, composed of equal and opposite halves. If P is what we 
call a complete reflector, the electric currents set up in it send waves 
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to the left as well as to the right and the nature of things is such that 
the waves which go to the right are out of step with the original set, 
and so cancel them, while the waves which go to the left cooperate with 
the original set to produce what we call reflection. If gravitational 
forces could be explained in terms of waves, and if matter were such 
that these waves, when impinging upon suitable matter, set up other 
gravitational waves, we might be able to realize a shield for gravitation 
in the same way that we shielded the radio, introduced at the beginning 
of this essay, against the waves from without by the metal container 
with which we surrounded it. 
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PONTRIAGIN’S MAXIMUM PRINCIPLE AND THE 
PRINCIPLE OF OPTIMALITY * 


BY 
CHARLES A. DESOER'! 


ABSTRACT 


This paper demonstrates that Pontriagin’s Maximum Principle may be de- 
rived from the principle of Optimality. It considers a control system described by 
x = f(x, u, 4) where the control vector u is restricted to a closed and bounded 
set. The optimal control steers the system from an initial state x) at to to 
a moving target in such a way that the cost of control along the optimal trajectory 

fol x°( 7), u°(7), is minimized. The optimal control satisfies the maximum 
principle at each point of the trajectory: max(§ (x), f (x, u)) = (Q(x), f (x, w°)) where 
Ue 


the vector g(x) is related to VI' where I'(Xo, x) is the minimal cost of going from Xp to x. 


I. INTRODUCTION 


Interest in optimal control problems has been growing at a very fast 
rate both in this country and in the Soviet Union (1-10).2. By optimal 
control we mean the control of a dynamical system such that the time 
integral of a function of the state and of the control signal(s) is mini- 
mized while the control signal is restricted to be in a specified closed 
and bounded region. As Pontriagin has shown (7), this formulation 
includes both the minimal time and the minimal fuel consumption 
problem. In 1959, Pontriagin (7) published a very interesting result 
known as the maximum principle. This result was particularly in- 
teresting because it stated properties of the optimal solution for very 
general nonlinear dynamical systems. Krasovskii obtained similar 
results (10). 

The purpose of this paper is to derive the maximum principle of 
Pontriagin in a completely different fashion. The point of view adopted 
here is to use the principle of optimality. As shown recently by S. E. 
Dreyfus (11), this point of view is particularly powerful to obtain the 
necessary conditions satisfied by the solutions of variational problems. 
As will become apparent in the following, the derivation is particularly 
simple and, in particular, does not require 2 to be convex. In fact our 
derivation will make it obvious that the validity of the maximum 
principle does not depend on the terminal conditions imposed on the 
optimal trajectory, that is, the maximum principle is valid if the opti- 


* The research reported in this paper has been sponsored by the Electronics Research 
Directorate of the Air Force Cambridge Research Center, Air Research and Development 
Command, under contract AF 19(604)-5466. 
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mal trajectory must reach a fixed point in state space, a moving point 
in state space (with a specified law of motion), or a moving subset of 
the state space. In their latest paper (8), Pontriagin et al. reached 
the same conclusion by extending their previous method. 


Il. STATEMENT OF THE PROBLEM 


We are considering a problem more general than that considered 
by Pontriagin (7). More precisely we consider a nonautonomous non- 
linear dynamic system S which can be controlled by 7 control signals 


Uy, U2, *-*, uy. The system is described by the law of motion of its 
state vector x in n-dimensional state space R": 
x = f(x, u, /) (1) 


where 


- x,) is the state vector in the n-dimensional 


(a) x = col(x1, - 
space R*; 

(6) the control vector u is restricted to remain in the closed 
and bounded set 2; for example, 2 might be the parallelepiped {u: | 1, | 
<yi,t =1,2,--- 

(c) the vector valued function f(x, u,¢) is assumed to have con- 
tinuous first partial derivatives at all points of R" X 2 X J, J being the 
real time axis. 


Any control function u(t) that satisfies condition (5) is said to be 
admissible. Suppose that x, the state vector of the system S, repre- 
sents the motion of an object in physical three-dimensional space. 
Let then x;, x2, x; be the Cartesian coordinates of its center of gravity 
in that space. The problem we are considering involves a target T 
represented by its law of motion in the physical three-dimensional space : 
z,(t), 22(t), 23(t), three specified functions of time. The problem is 
to steer the system S from the specified starting point Xo, to until it hits 
the target T (in physical three-space) while minimizing the integral 


where 


(1) fo is a specified cost function, 
(2) x(t) and u°(t) are the optimal trajectory and the optimal control, 
(3) ¢; is the instant of collision, that is, the instant at which x;(t,) 

= 2,(t:) fori = 1, 2,3. Observe that ¢, is not specified in advance. 


Before reformulating completely the problem in terms of concepts 
involving only the space R", observe that we may, without loss of 
generality, assume that f and fo, of (1) and (2), respectively, do not 
depend explicitly on time: indeed if they were dependent explicitly on 
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time we would introduce an (nm + 1) state variable x,,; = ¢ and an 
(n + 1)* differential equation z,,, = 1. Therefore, from now on we 
replace (1) and (2) by (1a) and (2a) respectively : 


x = f(x, u) xeR", ueQ (1a) 


The target motion in physical three space suggests the consideration 
of the cylinder C in R": C = {x(t)|x,(t) = 2,(t) for ¢ = 1, 2, 3; x; arbi- 
trary for 7 = 4,---,m}. In general, C is a subset of R" that moves 
with time. 

The problem can be restated in R” as follows: find an admissible 
control u(t), te[to, t:.) which will bring the system S from Xo, to on the 
cylinder C so that the integral (2a) is minimized. 


III. THE MAXIMUM PRINCIPLE 


Assume that an optimal trajectory exists. Let it be represented 
by x°(t), te[to, t:]. Consider an arbitrary point x°(¢) on the optimal 
trajectory, with to <¢ <¢,. For convenience we write henceforth 
x for x°(t). The principle of optimality states that an optimal control 
has the property that whatever the initial state and the initial decision 
are, the remaining decisions must constitute an optimal control with 
regard to the state resulting from the first decision. Let us express 
this fact for a special class of other possible trajectories. Let A be a 
small positive interval of time. Consider all admissible controls u over 
the interval (¢ — A,?t). From (1a), any trajectory using such a control 
will go through the point x — f(x, u)A + 0(A). If I'(Xo, x) is the mini- 
mal cost of going from X» to x (that is, the cost computed along the 
optimal trajectory), then the cost along the perturbed trajectory just 
mentioned is 


I'(Xo, x — f(x, u)A + 0(A)) + fo(x, u)A + 0(A) 


where the last two terms give the contribution to the cost given by 
(2a) during the interval (¢ — A,¢). Therefore the minimal cost of 


going from X» to x is characterized by 


T'(xo, x) = min {I'(xo, x — f(x, u)A + 0(A)] + fo(x, u)A} + 0(A) 


u(t)«2 


which expresses the optimality of the control u over the interval 
(t — A,t). Now if I'(xo, x) has continuous first derivatives in an open 
set containing x, we may use the mean value theorem for functions of 
several variables (12) and in the limit of A — 0 we get, after cancelling 
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terms and changing signs, 


0 = max { VI'(Xo, x) -f(x, u) — fo(x, u)} (3) 


uel 


where VI'(Xo, X) is the gradient of the scalar ['(xo, x) evaluated at x. 
In order to write condition (3) more elegantly, introduce the following 
n + 1 dimensional vectors 


f(x, u) = col (fo(x, wu), f:(x, u), ---, fa(x, (4) 


or oF 
0x1 Ox, 


(5) 


z(x) = col (- 1 


Hence, using the notation (a, b) to indicate the scalar product of a and b, 


max (2(x), f(x, u)) = 0. (6) 


uel 


This is the maximum principle. Its meaning is the following: if 
there is an optimal trajectory x°(t) (and an associated optimal control 
u’(t)) and if I'(xox)eC! at x, then at each point x of the trajectory the 
optimal control u° maximizes the scalar product in (6), that is, 


max (g(x), f(x, u)) = (g(x), f(x, u°)). (7) 


uel 


Furthermore, the maximum value of this scalar product is 0, 


(a(x), f(x, u°)) = 0. (8) 
IV. THE DIFFERENTIAL EQUATIONS SATISFIED BY 3(x) 


Consider the successive values taken by 2(x) along the optimal 
trajectory. We shall write 2(¢) for z[x°(t)]. For 7 = 1, 2, ---, n, we 
have, assuming that I has continuous second partial derivatives (12), 


r 
OX jOX; jai OX 


dgo 


since g) = — leverywhere. Now, observe that (8) defines u° implicitly 
for each x, that is, it provides a means for computing u® at each point x. 
Therefore as we consider (8) for various values of x, the corresponding 
value of u® varies with x and if we differentiate (8) with respect to 


UJ. 

, 

al 

and 
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x; we have 


OX; OX: OX; 


Hence the differential equations for the g;’s become 


at j=0 Ox; Ou, Ox; 


dgo 
dt 


Observe that if the maximizing u, namely, u° of (7) actually is in the 
interior of 2, then 


0 Zo(to) (9) 


n OF ; 
ou, 


and (9) simplifies to 
0 
j=0 
If, however, u° lies on the boundary of , this simplification will not 
occur. It is interesting to note that Eqs. 9 form the adjoint system 
to the first variation of (1a), provided in (1a) u® is considered as a 
function of x. 
Vv. REMARKS 

Equation 6 is a statement of Pontriagin’s maximum principle (7). 
Our derivation, however, is completely different and shows that the 
maximum principle holds irrespective of the terminal conditions im- 
posed on the optimal trajectory. 

The vector 2 satisfies also the adjoint system to the first variation 
of (1a), as does the vector J in Ref. 7. However, in the present deriva- 
tion u° appears as a function of x and, therefore, must contribute to the 
first variation. 

The above derivation could be repeated, with obvious modifications, 
if the minimal cost C(x, x:), of going from x to x; were used instead of 
I'(Xo, x). In general, VI! and VC bear no simple relation to each other 
because it is only along the optimal trajectory that (xo, x) + C(x, x:) 
= I'(X»,xX:)._ Examples show that in cases where (3) does not specify 
u® uniquely, the corresponding relation (with VI replaced by —VC) 
may determine u° uniquely. 

In the special case of a minimal time problem, that is, fo(x, u) = 1 
on R* X Q, the maximum principle reduces to 


max (Xo, x), f(x, u)) = (VI' (xo, x), f(x, = 


uel 
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Physically, this is obvious: the middle member of the above equality 
is the rate at which time increases along the optimal trajectory ; hence 
should always be equal to 1; the first member states that the optimal 
control u° should be such that at each point x, the optimal instantaneous 
velocity f(x, u®) should be directed so that the system reaches points 
furthest in time (as measured by I'(Xo, X)) as possible. 

In some cases the maximum principle points to significant properties 
of the optimal control. Consider, for example, a nonlinear autonomous 
system 

x = ¢(x) + Bu ueQ 
where 


(a) ¢(x) is a vector function of x, usually nonlinear, satisfying a 
Lipschitz condition for all x. 

(6) Bisann X r matrix. 

(c) Gis the parallelepiped |u;| < 4=1,2,---,7 

(d) fo(x, u) = 1 on R* Xr, that is, we consider a minimal time 
problem. Then the maximum principle states that at each point of the 
optimal trajectory, the optimal control maximizes 


(g, Bu) = (B’g, u) 
hence 


u; = sgn b,:2;) 


foi 


in other words, the optimal control for these nonlinear systems will be 
of the bang-bang type. 
VI. CONCLUSION 

We have shown that Pontriagin’s maximum principle simply ex- 
presses the optimality of the choice of control function at each point of 
the trajectory. The present derivation makes no assumption on the 
set 2 except that it is closed and bounded. In fact, the whole derivation 
shows that the maximum principle will always be satisfied irrespective 
of the particular terminal conditions imposed on the optimal trajectory. 

It should be kept in mind that the conditions derived here are 
necessary, that is, if an optimal solution exists and if f and [ have con- 
tinuous first derivatives, then the maximum principle holds at each 
point of the optimal trajectory. 
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A NOTE ON THE RELATIVISTIC BOLTZMANN EQUATION 
AND SOME APPLICATIONS 


BY 


S. YADAVALLI! 


4 


ABSTRACT 


Starting from the Boltzmann equation in the proper frame, the relativistic ver- 
sion of the Boltzmann equation in the absence of collisions is derived. As an applica- 
tion, a simple example of a velocity modulated drifting relativistic electron beam 
in the presence of a velocity distribution is treated. 


I. INTRODUCTION 


In the solution of many problems in plasma physics, the kinetic 
method involving the Boltzmann equation is often employed. The 
standard nonrelativistic form of the Boltzmann equation is well known 
(1).2. Relativistic versions of the Boltzmann equation have been either 
given or written in the past by others (2, 3, 4,5). An alternative deriva- 
tion of the Boltzmann equation in the absence of collision is given 
below, and the resulting equation is expressed in a form identical with 
that of Clemmow and Willson (2) for purposes of comparison. A simple 
example of a velocity modulated drifting electron beam is treated below. 


II. DERIVATION OF THE RELATIVISTIC BOLTZMANN EQUATION 


Let f(%, 3, t) be the one particle distribution function in phase space 
x, ) where X and @ are the ordinary spatial and velocity coordinates and 
t the time. Then, the Boltzmann equation in the absence of collisions 
may be written as 


df 
(1) 


Equation 1 is the standard non-relativistic form of the Boltzmann 


equation in the absence of collisions (1). 
In the relativistic case, one may start with the Boltzmann equation 


in the proper frame by writing 


af 
(2) 


! Stanford Research Institute, Menlo Park, Calif. 
* The boldface numbers in parentheses refer to the references appended to this paper. 
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and transform the equation to any other frame of interest, as done below. 
In Eq. 2, ds is the proper time interval times c (the speed of light), 


that is 


Equation 2 can also be written as 


af _ dx; Of du; _ 
ds ax; ds au; ds (3) 


where the summation convention of repeated indices is understood. 
In the four vector notation (6) 


x; = (x, y, 8, tc) = (x, y, 2, 7) 


= 


where Latin indices designate the four vectors, and Greek indices denote 


ordinary vectors; also 7 is the ordinary velocity and i = V—1. Note 
that: 7 = 1, 2, 3, and 4, and a = 1, 2, 3. One also has 


= 


d 


vy? \1/2 dt 
(1 


| 
; 
: 
dx ; 
3 
Ue = 
1 
1/2 
1 
du; 
2 1 
4 d 1 
y?\1/2 dt py? \1/2 
c(1-5 1-= 
2= 
dx,)? = — ds? 
= | 
= 
: 
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From the above relations, one notes that 


_ of Of for j= 1,2,3 


~ for j=4 
OX ; c 
(4) 


Substitution of relations (4) in the first term of Eq. 3 yields 


af dx; 1 af 5) 


V, denotes the spatial gradient. 
At this point we note the following. The distribution function f is 
a function of 7, v2, v3; f can also be expressed as a function of any three 
components of the four vector u; because conservation of energy and 
momentum are assured. Since it is desired to show the equivalence of 
this with the equation obtained in (2), we assume that the one particle 
distribution function is a function of 2, v2, and 03. 
From the above transformation relations, one has 


af due Ue 


1/2 2\1/2 
(4) 
c 


Substitution of relations (6) in the second term of Eq. 3 and addition 
of the same to (5) yields 


af d Ma 
@) 


where m» is the rest mass of the particle. 
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4 
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Note that no velocity dependent forces are assumed to exist in the 
derivation of Eq.7. It is clear that Eq. 7 of this paper, and Eq. 8 of (2) 
are identical but for a velocity dependent force term included by Clem- 
mow and Willson. It is not necessary to show the invariance of Eq. 7 
because the derivation presented here started from the proper frame 
and the Lorentz transformation is applied. However, if desired, a 
proof similar to that given in (2) may be stated. Noting that the ordi- 
nary force F is rate of change of momentum Eq. 7 may be written as 


(8) 


In the problems of interest, F is the Lorentz force and is given by 
e(E +% X B), where e is the magnitude of the charge, E and B being 
the electric and the magnetic field intensities, respectively. 

Note that in Eq. 8 the motion of positive ions is neglected. How- 
ever, to take account of the same one can introduce another distribution 
function g(%, 0, ¢) and write an equation similar to Eq. 8 above. If one 
assumes, however, that the positive ions are uniformly smeared out, 
one need not consider g(%, ¢). 

Although there are many applications of Eq. 8, we will consider 
below the problem of a drifting electron beam. One comes across such 
a problem in the discussion of microwave devices like the klystron. 


III. FORMULATION AND SOLUTION OF THE PROBLEM OF DRIFTING ELECTRON BEAM 


To simplify the mathematics only one dimensional electron streams 
are considered here. Imagine a parallel plane (diode type) region 
across which is impressed (or exists) a small alternating (a.c.) electric 
field. The distribution function then satisfies the Boltzmann equation 
in one dimensional form 


(1 
¥ 


af 
where ¢ is the electronic charge, x is the space coordinate normal to the 
planes, v is the x-directed velocity and E is the electric field. Let, 


E = | 
f fo(x, v) + filx, (10) 


| 
- 
v? \ 3/2 
; 
2 
+ ‘Vf = 
ot Mo 
where 
d 
2 1/2 
i—- — 
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so that f; «fo (fo the equilibrium distribution) and £; is assumed to 


be small. Substituting of relations (10) in Eq. 9 and linearizing (that 
is, neglecting the products of small quantities), one obtains Eq. 11 for 


the a.c. components of 
(1 y? 
a 


tof: + = (11) 
From Maxwell’s equations, one has 
dE, _ 


Integrating Eq. 11 one obtains 


filx, v) = f,(0, 


v 


(1 
+f v—E,(x’) ap dx’. (13) 


fi(x,v) can be determined from the integral equation (13) with a 
knowledge of v:(x), f:(0, v), fo(v) and Ei(x). The only boundary condi- 
tion needed is the value of f:(0,v). One can also obtain an integral 
equation for the a-c. current component as follows: 


in(x) = f° (14) 


In a one dimensional problem, with no sources inside the plasma (or 
electron stream), from Maxwell’s equations one obtains 


(15) 


1W€o 


Use of Eqs. 13, 14 and 15 yields the following integral equation for the 
a-c. component of the current, 


= ik vf (0, v) exp ( ) dv 


v 


x exp| (16) 


Vv 


7 ab 
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Drifting Electron Beam with Small-Amplitude Velocity-Modulation 


Although one can use the relativistic Maxwellian velocity distribu- 
tion (7) for fo(v) in the solution of the present problem we will not do so 
here. Instead for simplicity in mathematics which can be justified (8), 
a rectangular distribution for fo(v) will be assumed here. 

For the computation of f;(0, v), the model chosen is as follows. The 
stream is allowed to pass through two parallel grids across which there 
appears a voltage V, e*‘. It is assumed that no transit time effects 
occur and the grids are perfectly permeable, Vi « Vo, Vo being the 
voltage corresponding to the initial velocity of the electron stream. 

To the left of the gap, we have 


fo(v) = — — S(v —% — w)] (17) 
where po is the d-c. charge density at the entrance plane, w the width of 


the distribution function, and 


for v > V% 
S(v — ) = 0, otherwise 


(vo is the smallest velocity of an electron in the drift region defined by 


the d-c. potential.) 
Also, 


(+7) 


One can also write Eq. 17 in terms of voltages corresponding to the 
velocities as 


and 


folv) [S(V — Vo) — S(V — (17a) 


where V> and V,‘”) are the potentials (voltages) corresponding to 
velocities vp and (vp + w), respectively, and S is a step function already 
defined. 

To the right of the gap, we have 


fotfie = = {S(V — Vo — Viet) — S(V — Vo — Vie*)}. (18) 


Expanding the right hand side of Eq. 18 in a Taylor series and retaining 


eV 
only terms of first order in V; (7, = evr), one obtains 


: 

4 

: 

: 

V 
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v— — w) 


(vo + w) 


(19) 


— 
f:(0,v) Aa g(To) — g(ro™) ( 


where g(r) = [1 + r(1 + 7/2) ](1 + 7)- and 6 denotes a Dirac delta 
function. 

If it is further assumed that g(7ro‘”) ~ g(ro) (that is, a narrow dis- 
tribution is treated), then Eq. 19 further simplifies to Eq. 20 which will 
be employed in the ensuing calculations. 


— %) — % — w) 
”) Mo &(r0) Vo (vo + w) 


To obtain the a-c. component of current in a drifting electron beam 
which was initially velocity modulated, one has to solve Eq. 16 subject 
to the boundary condition given by Eq. 20. Equation 16 can be solved 
by Laplace transforms, and to do this introduce the Laplace transform 
defined by 


ti(p) = real part of p > 0. (21) 


Taking the Laplace transform of Eq. 16 and using the Faltung theorem, 
one obtains 


(1 


Molweo 


1 — 


It is obvious that 7;(x) can be obtained from 


1 fae, 
4;(x) = 4; (p)er*dp. 


We also have 


df o 


= — 09) — — w)}. 


Substitution of Eqs. 20 and 24 in Eq. 22 yields, after integration, 


4 

0 
p+ 
v 

0 
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where 


(26) 


3/2 2 
It is also assumed that (1 (1 ter) or, «1. 


0 
Finding the poles of 7:(p), pte neglecting terms of the order greater 
than one in w/o, after a little manipulation one finds that 


(1 - FS) 


x {sin (1 - exp| - =| (27) 


Vo Vo 


The non-relativistic analogue of Eq. 27 has been obtained previously 
in (8). 

Although it is already known, it is not out of place to draw the 
attention of the reader to two facts: 


(a) The plasma frequency that enters in is as defined by Eq. 26— 
where the so-called longitudinal mass of the electron is employed, and 

(b) The presence of g(r) (<1 in the relativistic case) indicates that 
it is much harder to modulate a relativistic electron beam than a non- 
relativistic one where g(r») is essentially unity. 


Finally it may be added that there are many other problems where 
the relativistic Boltzmann equation may be employed. 
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THEORY OF MAJORITY DECISION ELEMENTS 


BY 


SABURO MUROGA,' IWAO TODA AND SATORU TAKASU ? 


SUMMARY 


A decision-making organization whose output is either 1 or 0 according to which 
of the numbers, 1 or 0, is predominant in its inputs, is defined as the majority decision 
element. A Boolean function represented by such an element is called a majority 
decision function. Elements of this sort play a significant role in the computer field. 
The present paper deals with various problems concerned with a single element of this 
sort: that is, in Section 3, the algebraic properties, and a necessary condition for 
realizability of a given function; in Section 4, the types of functions realizable by a 
single element; in Section 5, the determination of the structure of the element for a 
given function by means of linear programming; and in Section 6, properties of the 
functions of a small number of variables realizable by a single element. 


1. INTRODUCTION 


The neurons in living organisms are believed to work on a majority 
principle and their behavior can be analyzed by propositional logic (1).* 
It has been shown that the majority decision principle reduces failures 
of neural signal transmission (2). 

Recently in electronic computers, logic elements based essentially 
on the majority decision principle, like magnetic cores (3) and para- 
metrons (4), have been used. Even circuits with vacuum tubes or 
transistors could be built on this principle and actually the principle 
had been used in the Kirchhoff adder. One of the authors has showed 
(5) that the idea of coupling numbers, that is, relative coupling ampli- 
tudes of inputs, makes even a single element based on the majority 
decision principle represent a fairly complex function and that a net- 
work for a given function can be realized with fewer elements of this sort 
and with a higher speed. As the switching speed of the physically 
realizable elements is limited, this realizability of a complex function by 
a single element may be of some importance for speeding up com- 
puters. Also interesting relations were shown between the probabilistic 
behavior of the element with its representable functions, and an error- 
reducing scheme of a network of such elements (6). 


11BM Research Center, Yorktown Heights, N. Y.; formerly with the Electrical Com- 
munication Laboratory, Nippon Telegraph and Telephone Public Corp., Tokyo, Japan. 

? The Electrical Communication Laboratory, Nippon Telegraph and Telephone Public 
Corp., Tokyo, Japan. 

3 The boldface numbers in parentheses refer to the references appended to this paper. 
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2. MATHEMATICAL REPRESENTATION OF THE MAJORITY DECISION ELEMENT 
Definition of the Majority Decision Element 


In the present paper, two states of firing and nonfiring of the neuron, 
and two different voltage values of a logic element in a computer are 
expressed by the binary values, 1 and 0. Of course there are other ex- 
pressions like +1 and —1. A decision-making organization whose 
output is either 1 or 0 according to which of the numbers, 1 or 0, is 
predominant in its inputs, will be henceforth called the majority deci- 
sion element. 

Generalizing this in this paper, a value weighted on the voting effect 
of each input value will be called a ‘‘coupling number’’ of the input. 
The majority decision will depend on this. Let coupling numbers 
of n variable binary inputs %1, %2, be Wi, We, Wn, 
respectively, and let w. be a coupling number of a constant input that 
takes permanently the value 1 or 0. Then an output value of the 
element is specified as 0 or 1 according to which of a sum of coupling 
numbers of inputs having 0’s and that of inputs having 1’s is larger. 
However it may be convenient to use the concept of a threshold T of 
the element as follows and to check whether a sum of coupling numbers 
of ‘‘variable”’ inputs having 1’s reaches this threshold value or not. 

w; (¢ = 1, 2, ---, m) in this paper is assumed to be a positive real 
number. When coupling of an input to an element is realized physically 
by wiring on a transformer as in the parametron, it is allowed only to 
be a positive or a negative integer (4). In a circuit with the Esaki 
diodes to which three-cycle excitation is applied (7, 8), coupling is 
realized with a resistor and it is a positive real number. However as a 
negative value of w,; can be reduced to a case of a positive w; by taking a 
negation of the input as a new variable, negativeness of w; is not 
essential as far as a single element of the majority decision principle is 
dealt with and will be henceforth eliminated. Thus, in this paper w; 
is limited to a positive real number and in most cases a positive integral 
value of w; will be considered. 


Definition 1: Let coupling numbers of variable inputs x1, 2, +++, Xn 
be wi, We, «++, W,, and a coupling number of a constant input be w.. 
The sum of these will be called the total coupling number of inputs: 


(1) 


Definition 2: A majority decision element is defined as an element 
whose output value is 


(2a) 


i=] 


< (2b) 


i=1 
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for n variable inputs x;, X2, Coupling numbers wy, we, ---, w, 
and 7 are assumed to exist and to be positive real numbers. 


When we can find a positive number e such that we, sT-e 


always holds as (2b), (2a) and (26) may be reduced to (4a) and (40), 
given later, by dividing (2a) and (20) by ¢ and replacing w;/e and 7//e 
by new w; and 7. Therefore (4a) and (4b) can be used henceforth in- 
stead of (2a) and (26). Note that the x,’s in (2a) and (2) are no longer 
variables but take fixed values of 1 or 0. 

The behavior of a majority decision element can be expressed by a 
Boolean function. For each subscript set a, 8, ---, y such that w. + ws 
+--+ +w, 2 T, construct a product xax,--- x, The logical sum of 
all of such products represents the Boolean function that expresses the 
majority decision element in Definition 2. This will be called the ma- 
jority decision function. Note that any product that includes another 
product as a factor can be deleted from the function, because the latter 
always takes the value of 1 when the former is 1. Thus we have a 
Boolean polynomial without negation of the variables. Henceforth 
throughout this paper we assume that a majority decision function ts in 
this irredundant form, that is, consists of only prime implicants. 

The structure of a majority decision element is completely specified 
by its threshold 7 and the coupling numbers of variable inputs, 
Wi, We, -**, Wa. The structure will henceforth be denoted simply as 
[W1,We, 7]. Alternatively, specification of w, instead of T also 
can specify the structure. It is easy to get a Boolean polynomial or the 
inequalities in Definition 2 from this structure, but in general, determi- 
nation of the structure from a given Boolean polynomial is not easy ; 
it is a major theme of the present paper. The threshold JT may be a 
number taken such that Sy < 7 < Sy. Here Sy is the maximum of 
sums of coupling numbers of the inputs, each of which does not reach 
the threshold with values of one of these inputs, and Sy is the minimum 
of those, each of which exceeds the threshold with values of one of 
these inputs. 

However, when all of w:, ws, ---, w, and w, are positive integers, 
and if the total coupling number of the inputs, K, is an odd integer, the 
threshold JT may be defined as a positive integer as follows: 


T = +.) (3) 


where — and + correspond to 1 and 0 of a constant input, respectively. 


Then > w,x, takes only positive integer values, so Eqs. 2a and 2b may 
i=1 


be replaced by the following: An output value of the element is 
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Lif Swe,-T20 


i=1 


i=1 
This paper will be based on this form of equations. Let concrete values 
of x,’s be expressed as x,;; and then the above may be written as follows: 
An output of the element is 
=1,2,---,7) (4a) 


i=l 


Relation Between the Physical Elements Based on Majority Dectsion 
Principle and Their Theoretical Models 


In the elements of Definition 2, variables take values either of 1 or 0 
and outputs are decided by the majority of these binary values. Physi- 
cally realized elements, however, like parametrons, neurons and Esaki- 
diode circuits, may not be interpreted to work concretely on values of 


1and0. It may be physically more accurate to state that parametrons 
work on the majority of +1 and —1 (4). In the case of neurons, it 
may be natural to regard 1 and 0 for excitatory inputs, and 0 and —1 
for inhibitory inputs, as the contribution of the majority decision (2). 
In the case of a neuron with both excitatory inputs and inhibitory 
inputs, an input may appear to take ternary values of 1, 0, and —1, but 
actually each of the inputs takes binary values, that is, 1 and 0, respec- 
tively for 1 and 0 of an excitatory input and —1 and 0, respectively, for 
1 and 0 of an inhibitory input. Then a theoretical model for the neuron 
where all inputs take the same binary values of 1 and 0 can be obtained 
as follows. When the number of inhibitory inputs is m, replace the 
threshold value T by T + m and take 0 and 1 instead of —1 and 0 of 
each inhibitory input, respectively, as the contribution to the majority 


decision. 
Situations similar to these appear in various circuits with the Esaki 


diodes. 

However, consideration of the binary values as being 1 and 0 is 
sufficient to treat theoretically the above physical devices. As shown 
in the following, it is easy to convert between values of input variables 
x,;’s in the theoretical model and those of the physical quantity &,’s 
that participate in majority decision of the output of the physical device. 

Suppose that £; assumes physical values either of a; or 6; correspond- 
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Then we have the following conversion formula 


=ax,+6,1 


ing to 1 or 0 of x,. 


(5) 


In the physical device based on the majority decision principle, its out- 
put assumes the value 


IV 


0, if w,é; T - 


Substitution of Eq. 5 gives the relation that an output of the element 
assumes the value 


1, if w; (a; b;)x; b, 
0, if w;(a; b,)x; < b; — 1. 


This represents a theoretical model with a threshold 7 — > 6; and with 
coupling numbers w;(a; — 0,;) for x;’s, corresponding to the above de- 
vice. Conversely, the structure of the physical device can be specified 
from that of the theoretical model if a; and 6; of the device are known. 

Consequently, the theoretical model of Definition 2 is sufficient 
for discussions of the logical behavior of any physical devices based on 
the majority decision principle. 


3. SOME PROPERTIES OF THE MAJORITY DECISION FUNCTION AND A NECESSARY 
CONDITION FOR ITS REALIZABILITY BY A SINGLE ELEMENT 


Real numbers named as coupling numbers are attached to inputs of 
the majority decision element and ordinal relations between them 
specify the behavior of the element. Here we shall discuss how a rela- 
tion between coupling numbers and a threshold is reflected in a property 
of the majority decision function itself. 

An inclusion relationship between Boc'ean functions may be de- 
fined as follows: Consider two functions f(x,, ---, X,) and g(x1, «++, Xn); 
if any set of x1, x2, ---, x, satisfying f(x1, ---, x.) = 1 satisfies also 
g(x1, ---,X,) = 1, but the converse does not necessarily hold, we write 


Vids, ***, Xe) (6) 
In addition, if Eq. 6 holds and, further, there exists a set of x1, ---, x, 
satisfying simultaneously f(x,, ---,x,) = 1 and g(x, ---,*,) = 0, we 
write 


A function f(x, ---,x,), where for any 7(¢ = 1, 2, ---, m), 
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holds, is called a positive function. A positive function can be ex- 
pressed by a polynomial of xi, ---, x, without negations of these vari- 
ables and vice versa. A majority decision function is a positive func- 
tion but the converse is not true (5). The positive functions defined 
here are the same as E. N. Gilbert’s switching functions that are ex- 
pressed without negations of variables (9). 


Ordinal Relations Among Input Variables 


A structure of the majority decision element specifies the function 
that the element represents. However, the structure of an element 
that is to represent a given majority decision function is not unique. 
Consider the following example. The majority decision function 
x; + x2 is represented by the element with the structure w, = w. = 1 
and T = 1 and also by w; = 2, w. = 3 and T = 2. However, it may 
be important to see as in Theorem 2 that when a majority decision 
function is given, at least ordinal relations among coupling numbers 
can be uniquely determined with respect to any one of the elements 
that are to represent the given function. The following definition will 
be useful. 


Definition 3: Consider a positive function f(x1, 2, -:+,%n); a re- 
lation 

ij 
ten 2 f(*1, 


f (x1, (9) 


for any set of 7 and j( # 2), is denoted as 


xX (10) 
Verbally it may be stated as “‘the input variable x; is not smaller than 
x;.’ In particular, when an equality is not included in (9), it is denoted 
as x; > x; and when only equality is included there, x; ~ x;. For com- 
pleteness of definition, x; ~ x; is defined. 

This binary relation has the following properties: 


> on the positive functions is a par- 


Theorem 1: The binary relation 2 
tial ordering. 


Proof: Obviously this binary relation satisfies the reflexive and 
anti-symmetric laws from the definition. Assume x; 2 x2 and x2 2 Xs, 
that is, 


(11 


; 
Bt 
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From a property of the positive functions 


=f(1, 1, 0, x4, + f(1, 0, 0, x4, 
+++) = f(0, 1, ++ + f(0, 0, 1, ° 


hold. Inserting x; = 1 and x; = 1 in the two inequalities of (11), 
respectively, we have f(1, 1, 0, x4, --+, x.) Df(1, 0, 1, x4, +++, Xn) 
> f(0, 1, 1, x4, --+, x,), and thus the first term in the first expansion in 
(12) is not smaller than that in the second expansion. Similarly, the 
second terms are compared by inserting x; = 0 and x; = 0 in (11). 
Consequently, as the terms in the first expansion of (12) were shown to 
be not smaller than the corresponding terms in the second expansion, 
we have f(1, x2, 0, x4, DF(O, x2, 1, x4, «++, Xn), that is, x1 > x; 
and the transitive law is satisfied. 

When we proceed into the majority decision functions from the 
positive functions, this relation will prove to be a total ordering. 


X2, 0, X4, 
F(O, x2, 1, X4, 


(12) 


‘) 
) 


Theorem 2: The binary relation = for a majority decision function 
(x1, X2, +++, X,) is a total ordering. w,’s that have an ordinal relation 
among them as real numbers are homomorphic with the input variables 
x,’s defined with the binary relations >. That is, if w; > w,;, then 


(13) 


and if w; = wy;, then 
(14) 


Proof: As the relations were shown to be a partial ordering in 
Theorem 1, if either x; 2 x; or x; < x; for any set of x; and x; holds, 
the relations are total orderings. If the majority decision function 
(i, X2, +++, Xn) isrealized with a single element of a structure [w1, we, ---, 
w,; Tj, one of the following three relations must hold between coupling 
numbers w; and w; of any two variables: 


(i) When w; > w,;, if for a set of values of variables excluding x; 
and x; 


> wx 2 T — w,, (15) 


t j 
holds, f(*1, 1, ---,0, +--+, %,) takes value 1. And if 
wx, 2 T (16) 


1 j 
holds, f(x;, ---,0,---,1, +++, %,) takes value 1. Since T — w; 


7) | | 
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< T — w; holds from the assumption, a set of variableZvalues 
(excluding x; and x;) which satisfies (16) satisfies (15) also and 
therefore 


f (x1, -++,0, Xn) > f (x1, »++,0, Xn) (17) 


holds. Consequently we have x; 2 xj. 


(ii) When w; = w;, we can conclude by a discussion similar to that of 
case (i) that only an equality holds in (17) and therefore x; ~ x;. 


(iii) When w; < w,, we again, as with case (i), have x; < x;. 
Consequently, since for any x; and x; a binary relation 

< exists, the homomorphism of x; defined in this way with w,’s 

as real numbers has been proved. 


From this theorem, input variables x,’s can be arranged in order by 
means of the binary relation 2. Even though the structure that is ac- 
tually required to realize a given function is not known, ordinal relations 
among coupling numbers to be realized can be seen from the binary 
relations peculiar to the function. Of course if x; > x;, then w; > w,; 
holds. For example, for the function f(x1, x2, = %1 + as 
f(1, 0, xs) = 1 DF(O, 1, x2) = x3 and f(x1, 0,1) = x1 = f(x, 1, 0) hold, 
we have x; > x, ~ x3. This function is actually realized with an ele- 
ment of structure [2, 1, 1; 2] where there are the relations w; > wz = Ws. 

The following corollary is convenient to check whether a given posi- 
tive function is actually a majority decision function (17). 


Corollary 1: A majority decision function f(x:, ---,%,) can be 
written in the following form: 


t 


f (x1, Xn) = M (x:, es 


where x; is any one of the largest variables compared in the binary 
relation >. Here a relation 


MDN (19) 


1 
holds. Both M(x, -- and N(x1, Xn) represent 
the majority decision functions of (nm — 1) variables excluding x;. 
These can be further expanded as follows: 


i 
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j 
Xn) X; 


Xn), 


with respect to a certain x;(j # 7) that exists always and the relations 


(21) 


Mm, ) Ne 


j 
hold. , °**,Xn) etc. mean functions of (n — 2) 
variables excluding x; and x;. 


If we assume w; > w; > ---, the proof is easily obtained since f 
consists of prime implicants. 


Binary Relations among Input Vectors 


In the previous section the binary relations are defined by giving 
concrete values of 1 and 0 to only two of n variables. Here, generalizing 
this, 1 and 0 will be given to all m variables. While in the previous 
section coupling numbers of input variables are compared pairwise, the 
treatment in this section implies comparison of sums of coupling num- 
bers of input variables having value 1. When binary values are con- 
cretely assigned to all » variables, those sums of coupling numbers of 
variables having value 1, which exceed the threshold, specify the 
majority decision functions for the element. Though the ordinal rela- 
tions among these sums, for various sets of values of variables, can be 
easily determined if a structure of the element is known, they can be 
specified to some extent from properties of the given function even if it 
is not known, as shown in Section 3 of Theorem 3. However as shown 
in Section 2 of Theorem 3, this relation is not a true ordering, but a 
“pseudo ordering.” 

Concrete values of input variables x,’s will be assigned by &,’s and a 
particular set of these values by a vector expression = = (&, ---, &,). 
This will be called an ‘input vector.”” The j* vector in these 2" input 
vectors will be written as: 


where j = 0,1, ---, 2" — 1 and = Oor1 = 1,2, -+-, m). 


: 

Jj 
~ ~ 
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Definition 4: Let =; and Z: be two input vectors, where m out of n 
components of £;“) and &;® (¢ = 1, 2, ---, m) are of the same values. 
For example, 


§,@) =m+1,---,n) (23) 


(in general z’s for these equalities and inequalities are not necessarily 
consecutive). If 


holds for two functions of x,, ---, X» out of the m variables, this fact will 
be denoted by 


(25) 


where equality and inequality in (24) and (25) correspond to each other. 
Though this binary relation may be reduced to = (see Definition 3) 
f 


if only two components differ in £; and £;®, a different notation > is 
used here to stress a generalized idea. Such variables as x;, ---, X» in 
this definition may be called free variables and the rest constrained vari- 
ables. For example, when f(x;, %2, = + + %3%1 is given, 
for Z, = (0, 1, 0) and =, = (1, 1, 1), x2 is a free variable and x; and x; 
constrained variables. As f(0, x2, 0) = 0 and f(1, x2, 1) = 1, we have 


f 
1< 2s. (26) 


Theorem 3: If f(x:, -++,%,) is a majority decision function, the 
following holds: 

1. A binary relation is always determinable between any two of the 
2" input vectors. 

f f f 

2. If Zi < Ze and Z. < Zs, then Z, < =; (note that this is not a 
truly ordinal relation). 

3. In the special case where a structure [wi, ws, ---, Wa; J] of a 
majority decision element realizing this function is known: if } w,& 


> ¥ then 
i=1 


f 
12> (27a) 


and if wt, = then 


i=1 


‘ 
Ei > Es, = 
4 
: 
= 
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Proof: Let us start from Section 3 above. In the comparison of 
=, and Zo, assume x, ---, x, to be free variables and the rest to be con- 


strained variables. As a structure of the element i is given, 3 w+: for 


j = 1 and 2 can be calculated. Since = w, holds, if 
i=1 

one of >, = and < exists between the remaining > wyé, and 
n i=m+1 

X< wt, then one of D>, = and C holds, respectively, between 


i=m+1 

Fity, +> +; Sea texi™, and f(x1, +++, Emai®, &,@), by 
a treatment similar to the proof of Theorem 2. Thus Section 3 of 
Theorem 3 is proved. Next, if it is a majority decision function, there 
should exist coupling numbers of an element to realize it, for which the 
property described above holds for any set of =, and Z.. Then always 


f 

a binary relation > can be established and Section 1 of Theorem 3 is 

proved. Last, for a structure of an element to realize ff (x1, °° 11% ), if 
f 

=1 < Ze and 2:2 =3 hold, we have W; Mi< x < wt, 


i=1 
from Section 3 of Theorem 3. Thus 2, <x =; is eS proving Sec- 
tion 2. 

From this theorem, the structure of an element realizing a given 
majority decision function can be seen to some extent. For example, 
assume that f(x1, ---, = + is given and compare 

= (1,1,0,0) and =. = (0,1, 1,1), where only is a free variable 
and the rest constrained variables. Then f(1, x2, 0,0) =x. and 


f 
f(0, x2, 1,1) = 0. As x. D0, (1, 1, 0,0) > (0,1, 1, 1) holds. From 
Section 3 of Theorem 3, if this function can be realized with a major- 
ity decision element, a relation w,; > ws; + ws must hold among its 
coupling numbers. Actually a structure of an element for this function 
is (3, 2, 1, 1; 5] (16). 

When there exists only one constrained variable, this theorem is 
merely a definition of the positive functions. When Theorems 2 and 3 
are compared, their physical implications may be different, as the former 
compares individual coupling numbers pairwise but the latter compares 
various sums of coupling numbers. However, if regarded as a neces- 
sary condition for realizability of a given function by a single majority 
decision element, the latter is a stronger statement because the former 
always holds if the latter holds, but the converse does not necessarily 
hold. In Theorem 2 only a binary relation between any two variables 
is a necessary condition and its satisfaction of the ordering axiom is a 
necessary condition. Let us give an example of a function that does 
not satisfy Theorem 3 though it satisfies the condition of Theorem 2. 
For f (x1, Xs) = { (x3 + + + (x4 X6)X3 + 
+ + + relations x; > x, ~ ~ X4 > Xs ~ Can 


5 
2 
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be easily ascertained. But as f(1, 0, 0, x4, Xs, %6) = XaX5 + X5X— + Xero 
and f(0, 1, 1, x4, %s, Xs) = x4 hold, there does not exist an inclusion rela- 
tionship between these two. Thus this function satisfies Theorem 2 
but does not satisfy Theorem 3. Also a necessary condition obtained 
previously (5) compares actually individual coupling numbers and then 
is similar to Theorem 2. As shown later in the next section, there 
always exists an inclusion relation between a function satisfying Section 
1 of Theorem 3 and its dual function. Finally, a necessary condition as 
the strongest among the ones obtained so far may be stated as follows: 
A necessary condition for realizability of a given function by a single 
majority decision element is Sections 1 and 2 of Theorem 3. 


Duality of the Majority Decision Functions 


When a function f(x:, ---,%,) is realized by an element with the 
structure ---, 7], its dual function f* (x1, ---, x,) is realized by 


an element with the structure [w, ---, w.; 1 — T + X w,], by negat- 
i=1 


ing a value of a constant input (that is, if it is 1, convert it into 0 and 
vice versa) without change in coupling numbers of the constant input 
and variables (5). The element represents the same function since the 
element works on the majority of a number of 0’s instead of that of 1’s 
when all inputs and an output of the element are negated, that is, 1’s 
and 0’s are reversed. Now we will show that in any function satisfying 
condition 1 of Theorem 3 and further in the majority decision func- 
tions there always exists an inclusion relation between f and f*. 


Theorem 4: Consider a function f(x., ---, x,); if the binary relation 
defined in Definition 4 holds between any two of 2* input vectors, 
either f D> f* or f C f* holds. 

Proof: A vector where all components of an input vector = are 
negated will be denoted by Z (forexample, Z = (0,1, 1) for Z = (1,0, 0)). 
Let A be the set of all input vectors which gives f the value of 1, and B 
be the set of all input vectors which gives f the value of 1, where f means 
negation of f. Also, let A* and B* be similar sets of input vectors for f* 
and f*, respectively. Define A’ = {Z|ZeA} and B’ = {2|ZeB} and 
then from the definition of the dual functions we have 


A* = B’ and B* = A’. (28) 

Now we will show that both 
Af\A' #¢ 


do not hold simultaneously if the assumption of the theorem is satisfied, 
where ¢ means the null set. Assume that both hold simultaneously. 


(29) 


: 
New 
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Then there should exist input vectors Z, and Z, such that Z,eA and 
=,eA, =.<B and =.eB. As the free variables are obviously the same in 
the comparison of =, and =, as in the comparison of Z, and Z:, decom- 
position of the components in these four vectors leads to =; = (X, ®), 
=, = (X, 6), =. = (X, Vv) and =, = (X, ¥), where X is a vector con- 
sisting of the free variables, # and W are vectors consisting of the con- 
strained variables and X, 4 and © are their component-wise negated 
vectors. As all components corresponding to the constrained variables 
in =, and &, differ, we have € = ¥. From the assumption, binary 
relations hold between =, and =, and between =, and 2», and for vectors 
belonging to A and B we have always f = 1 and f = 0, respectively. 
Consequently, with X representing agreeing variables, there must hold 


f(X, DF(X, (30) 
f(X, 8) Df(X, 


without equalities. However, as this contradicts the relation # = ¥, 
at least one of assumed relations must be an equality. If 


BNB'=¢ (31) 
is assumed, A and B divide 2" input vectors into two parts and then we 
have 

(32) 
With (28) this becomes A D> A*. Then f Df* holds. Similarly, 
f Cf* for A 1) A’ = ¢. 
For the self-dual functions, A (1) A’ = @ and B Qf) B’ = ¢ in (29) 
hold simultaneously. 
From this we have the following observation. 


Corollary 2: Between a majority decision function f andi ts dual 
function f*, either f > f* or f C f* always holds. 


Variation of the Threshold 


Here we shall discuss how the majority decision function repre- 
sented by an element is affected by variation of its threshold, while 
fixing the coupling numbers of variable inputs. For simplicity, in this 
section all coupling numbers and the threshold are assumed as integers. 
A majority decision function for an element with a threshold t¢ will be 


denoted as f;(x1, Xn). 
Theorem 5: lf = j, fi(x1, +++, Xn) S fi (x1, «++, Xn) holds. 
Proof: Only when x,, ---,x, satisfy 2 f; takes the 
k=l 


value of 1. Then from the assumption, 5 wx, 21 = 7 is satisfied 
k=1 


and for these f; = 1 holds. 


ie 
; 
4 
: 
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; 
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Theorem 6: The number of different functions that are realized by 
changing the threshold of an element but fixing the coupling numbers 
of variable inputs, is not greater than 2*-+ 1, including two functions 
which are permanently 1 and 0. 


Proof: There are 2” — 1 different sums of coupling numbers, 
Wi, We, Wap Wi + We, Wy +w,. These sums 
can be arranged according to their magnitudes and in general some of 
them may be equal. This greatest number of different sums can be 
actually obtained as follows. If we set w; = 2‘—', then the successive 
magnitudes each differ by 1 and then variation of the threshold within 
this range gives at most 2" — 1 different functions because a different T 
gives a different bisection of 2" — 1 sums, resulting in a different func- 
tion. Here we have the function 0 for 7 = 2* and the function 1 for 
T = 0. In the case of w; = 2‘-', the forms of functions for a given T 
can be easily obtained. 

For example, set w; = 1, w, = 2, and w; = 4. If we increase J one 
by one, we have functions x; + x2 + %3, X2 + XiX2 + Xs, XiXs 
+ and accordingly. (Note that only the function 
X1X_ + X1X3 + Xex; is excluded from all the majority decision functions 
of three variables.) However, in general any one of these functions can 
be realized with smaller coupling numbers. 


Convolution of Structures of Two Majority Decision Elements 


In this paper, convolution of two elements means making a majority 
decision element with coupling numbers w; + wy’, ---, W, + w,’, from 
the one element with wi, ---, w, and the other element with wy’, ---, w,’. 
In this section all coupling numbers and thresholds are assumed as 
integers. 


Theorem 7: Let mi(x1, +++, Xn) and m.(x1, X,) be majority deci- 
sion functions of two elements with two sets of coupling numbers, 
Wi, ***,W, and wy,’, ---, w,’, respectively, and threshold ¢. The ma- 
jority decision function, represented by an element with a threshold T 
which is made by convolution of these two elements, can be written as 
follows: 


where M = > w, and where, for t s 0, m, and n, are regarded as the 
i=1 


function which takes the value of 1 permanently. 


Proof: The polynomial fr is a sum of products, each of which 
consists of variables with values of 1 in an input vector such that 


4 
‘ 
t=0 
‘ 
: 
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> (w; + w,’)x; 2 T is satisfied. Let S be a set of such input vectors 


and a sum of w,’s of components for any of these vectors be ¥ w,x; = t, 
i=l 


where ¢ is within the range 0 <¢ S M. For the remaining w,’’s, 


>< wx; 2 T —t must hold. Then m-_,-n, takes the value of 1 as 


both mr_, and n, are 1. Conversely, only when both ¥ w,x; 2 ¢t and 


i=1 


x wax; 2 T —1t are satisfied, mr_,-n, takes the value of 1 and all 
i=l 


input vectors validating these inequalities belong to S. Consequently, 
M 


frand > mr_,-n, are simultaneously 1 for input vectors in S and simul- 


taneously 0 for input vectors outside S. Thus expansion (33) is 
obtained. 

An inductive synthesis of certain types of functions with n variables 
from a function of (m — 1) variables realized by a single element may 
be obtained as follows. 


Corollary 3: If a function f,(x1, «++, Xn-1) of (m — 1) variables real- 
ized by an element with fixed coupling numbers is known for any 
possible threshold value ¢, a function of m variables 


fu(x, 


can be realized with a single element, where between arbitrary thresholds 
u and v, a relation u > v holds. 


(34) 


Proof: Set (@=1,---,n —1), w, #0 and w,’ #0 

(¢ = 1,---,m —1); w,’ = 0 in Theorem 7. Then we have my = 1 

and for ¢ such that M =w, 2t 21, m, =x,. Thus, rewriting nr 
M 

and mr_, in (33) as fr and fr_., we have fr(x1, Xn-1) + 
t=1 


(x1, As fr_w > for t < M holds from Theorem 5, this 
equation is simply 


fr + (35) 


Though w,’’s (¢ = 1, ---, m — 1) are given, values of w, and JT may 
be optionally selected and therefore the stated function is realizable. 
As a special case of the majority decision functions of (34), set 
T = M and then +x, is realizable. If we set 
n—-1 
T = Dw,’ + 1, we have x,f(x1, «++, Xn-1) as mr = 0 holds. Also for 


i=1 


t 
n 
: 
n 
; 
4 ; 
4 
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T= M+1andT = & wy’, f(xs, + (41 + + n 


and x1 + +++, Xn—-1), respectively, are realizable. 


4. FORMS OF FUNCTIONS REALIZABLE BY A SINGLE MAJORITY 
DECISION ELEMENT AND VARIOUS BOUNDS 

It is an interesting problem to find how many functions of n vari- 

ables can be realized by a single element, by adjusting coupling numbers 

and a threshold. Here upper and lower bounds are discussed since 

exact enumeration seems difficult. 


Forms of Realizable Functions 


Though enumeration of realizable types of functions by a single 
element gives its lower bound, simple types of functions will be dis- 
cussed in this section. As shown already (5), 


Theorem 8: A symmetric polynomial of degree d for n variables 
can be realized by a single majority decision element of structure 
Wi =We=--- =w, = 1andT = d. 


The following function, where certain terms are eliminated from the 
above polynomial, is also realizable. 


Theorem 9: A function in which a symmetric polynomial of the 
second degree for m(< n) variables is eliminated from a symmetric 
polynomial of the second degree for m variables, can be realized by 
a single majority decision element. The number of such functions 

is (7 ), where —-12m23. (A func- 
m =3 

tion for m = 2 is a special case for the next theorem with d = 2.) 


Theorem 10: A function where one term is eliminated from a sym- 
metric polynomial of degree d for nm variables can be realized by a 
single majority decision element. The number of such functions 


Proof: Theorems 9 and 10 will be proved. Eliminate a sym- 
metric polynomial for m out of n variables from a symmetric poly- 
nomial of degree d for n variables. We get a realization as follows: 
Let a be the coupling numbers of the variables contained in the elimi- 
nated polynomial and (« + 8) be the coupling numbers of the remaining 
(n — m) variables, where a and 6 are positive integers, and let the 
threshold of the element be JT = am +1. Then a structure of the 
element is specified. First, since a sum of coupling numbers of m 
variables included in the eliminated polynomial is at most am and does 


4 
: 
: 
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not reach its threshold, a product of these m variables is not included 
in the function for the element. Next check a product of (d — k) 
variables out of these m variables and k variables out of the remaining 
(n — m) variables, where naturally n >m 2d 2k 21. A sum of 
coupling numbers of variables that are included in this product of 
degree d is 


(a + B)k + a(d — k) = ad + Bk. (36) 


If 


B =a(m—d)+1 (37) 
is chosen, Eq. 36 becomes a{(m — d)k +d} +k = T= am+1 as 
m — d = Oand the sides are equal for k = 1, and then this value reaches 
its threshold. It remains to be shown that when any one variable is 
eliminated further from such a term, the sum of coupling numbers of 
the remaining variables does not reach the threshold. The largest sum 
of coupling numbers of variables of this term where one variable is 
excluded is 


(a + B)(d — 1). (38) 


In order that this be smaller than 7, 


(d —1){a(m —d +1) +1} Sam 


(39) 


Then we have 


must hold. 


(40) 


This will be satisfied by setting 
a=d-1 (41) 


for m = d and particularly for d = 2 by setting a = 1, irrelevant of m. 
Thus a + 8 and a are d and d — 1, respectively, in the former case and 
m and 1 in the latter. The above proof is not valid for d = 1 but it is 
easy to see that Theorem 10 holds even for this case because the func- 
tion is of the first degree. 


Theorem 11: A function obtained from a symmetric polynomial of 
the second degree for ” variables by eliminating a symmetric poly- 
nomial of the second degree for m variables and adding a symmetric 
polynomial of not less than the third degree for the same m variables, 
can be realized by a single majority decision element and this number 
is given by 


= ) im (42) 


m=3 


where n — 1 = 3. 


ai 
es 
; 
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Theorem 12: A function in which a symmetric polynomial of the d‘* 
degree for m variables contained in a symmetric polynomial of the 
d* degree for n variables is replaced by a symmetric polynomial of the 
(d + 1)* degree for the same m variables, can be realized by a single 
majority decision element. This number is given by 


= 


wheren 


Proof: Both theorems above will be proved together. a will be 
the coupling numbers of m variables and a + 6 the coupling numbers of 
(n — m) variables and the threshold of the element will be T = ad 
+ af + 1, where a and £ are positive integers and f a non-negative in- 
teger. With respect to the first m variables, any product of (d + f + 1) 
variables of these has a sum of coupling numbers equal to a(d + f + 1) 
which is over the threshold, but any product of (d + f) variables has a 
sum of coupling numbers which does not reach the threshold. There- 
fore a symmetric polynomial of (d +,f + 1)'® degree for m variables 
exists in the polynomial which the element represents. Next, a sum of 
coupling numbers corresponding to a product which consists of k of the 
remaining (n — m) variables and (d — k) of the first m variables, is 
ad + Bk, being calculated in a manner similar to Eq. 36; it always ex- 
ceeds the threshold for k 2 1 if we choose 


B = of +1. (44) 


The largest sum of coupling numbers for a product in which a single 
variable is eliminated from this product is the value equal to (38). 
In order not to reach the threshold 


(45) 


(a+ B)@ — 1) Sad +af 


must hold. Substitution of Eq. 44 gives 


alf(d.+ 2) (46) 


Then when f = O ord = 2 holds, (45) holds if we seta = d-—1. Thus 
the theorems are proved. a+ Sand aaref + 2 and 1, respectively, in 
Theorem 11 and d and d — 1, respectively, in Theorem 12. It is easy 
to show similarly that the latter theorem holds even for d = 1. 

An example will be shown. The symmetric polynomial of the second 
degree for four variables 


+ + + + + (47) 


it 
: 
: 
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can be realized by a single element. According to Theorem 9, make 
the following two functions by eliminating a symmetric polynomial of 
the second degree for m = 2, x3x,and one for m = 3, xox + Xox4 + Xs%u, 
respectively, from this function: 


and 
XyXo + (49) 


Then both can be realized by a single element. Also, according to 
Theorem 10, any function made by eliminating any single term from 
(47) can be realized by a single element. According to Theorem 11 or 
12, the function 


in which a symmetric polynomial of the second degree for m = 3 in 
(47) is replaced by the third degree product, x2%3x,4, can be realized by 
a single element. 


Theorem 13: A function in which a single variable in a symmetric 
polynomial of the d* degree for n variables is replaced throughout by a 
product of any number of new variables can be realized by a single 
majority decision element. 


Proof: The variable x, is to be replaced by a product of m vari- 
ables Ym. Multiply all the coupling numbers of X,-1 and 
T by m. Coupling numbers of y:, ---,and y, should be set to one. 
Realization of the function is then obvious. 

For example, if we replace x3 in x:%_2 + %1%3 + X2%3 by xs%4, we have 
a realizable function 


In an asymmetric polynomial, however, such replacement does not 
always give us realizable functions. 
Lower Bound on a Number of the Majority Decision Functions 


In this and the following sections, a number of the functions that 
can be realized by a single majority decision element will be discussed. 
First the symmetric polynomials in Theorem 8 will be considered. 


There are ( b combinations of k out of m variables and the degrees 


of symmetric polynomials of k variables range from the first to the k**. 
Negation of different sets of variables gives different functions, Thus 
the number of such functions is 


: 
3 
: 
: 
4 
= 
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2-k-(7)| +2, (51) 


k=l 


including two functions that take constant values of 1 and 0. Summing 
the contents of the brackets, we have the following theorem. 


Theorem 14: The number of symmetric polynomials for at most m 
variables that can be realized by a single majority decision element is 


+- 2, (52) 


including functions made by permutations and negations of variables. 

This number is very small compared with 22”, the number of func- 
tions of m variables. However, we shall show that the number of asym- 
metric polynomials is much greater, though it is still much smaller than 
22" for large n. Functions shown in Theorems 9 to 13 are asymmetric 
polynomials based on symmetric ones but their number is quite small. 
The structures of two types of asymmetric polynomials in the next 
theorem are entirely different from those above, but their number is 
large enough to deduce a reasonable lower bound on the number of 
functions. Forms of these functions in general are complex and it may 
be difficult to show them explicitly, but their differentiation is possible. 


Theorem 15: A lower bound on the number of asymmetric poly- 
nomials for at most » variables that a single majority decision element 
can represent, is 


(26 oom) + 


n log n 
+29 (Cn/log n] — (53) 


for n = 9, where [x] means the integer that is the nearest to but not 
smaller than x (this notation is different from the one used in mathe- 
matical text books) and where log m is the natural logarithm of n. 


Proof: We shall discuss first the type of the functions charac- 
terized by having only a single product of 7** degree and other products 
of degree lower than 7. Henceforth the threshold T of the element will 
be kept constant, by adjusting properly a coupling number and a value 
of the constant input. Now, keeping coupling numbers of T of n 
variables, for example, x1, x2, «++, xr, as 1, independently vary the 
coupling numbers of the remaining (” — 7) variables through integral 
values from 2 to T. Th: function represented by this element always 
consists of a product of 7** degree x:x, --- x7 and products of lower 
degree, each of which includes at least one of the other variables. 


= 
: 
: 
: 
: 
: 
i 
: 
4 
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Since the 7“ degree product term cannot be deleted from the function 
expression by including any other term as its factor, all functions for 
this 7 differ from those obtained by starting from different values of 
T. Let w, be a coupling number of some variable x, of the remaining 
(n — T) variables, where 7 2w, 22. A product of this x, and 
(T — w,) variables of x;, ---, x7 exists in this polynomial function and 
a sum of their coupling numbers takes a value of just 7 when all factor 
variables take the value of 1. For a different value of w,, the number 
of variables in this product extracted from x, ---, x7 differs. Now 
if we compare two functions for any two different sets of coupling num- 
bers for variables x7,, to x,, we can find at least one different term such 
as the above. Since all possibilities produce different functions, the 
number of polynomials made in this way is 


(T — 1)*-*. (54) 


Selection of T variables and variation of T from 2 to n — 1 give poly- 
nomials of a number 
n 


(T — 1)*-?. (55) 


Here the case T = n is omitted since it represents a symmetric poly- 
nomial. Consequently the number of such polynomials for variables 
of a number equal to or less than which can be represented by a single 


element is 
> b (T — 1)*-T], (56) 
k=l T=2 


taking account of negation and permutation of the variables. Though 
evaluation of this formula is difficult, (7 — 1)*~7 has a larger effect on 


k 
the magnitude of the formula than ( fy for large kand 7. Examina- 


tion of the largest effect of this leads to the conclusion that by choosing 
one term of the double sum: 


T = [n/log. n ], k =n, (57) 
formula 56 is bounded below by 


n —[n/log n 


where [x ] is an integer which is the nearest to but not less than x. 
Next we shall discuss the second type of functions. These are ob- 

tained by changing the threshold value after assigning consecutive 

integers from 1 to m to coupling numbers of m variables as follows: 


4 
3 
: 
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w, =1,---,w; =4,+--,w, =m. There are some cases where we 
can have explicit forms of these polynomials for some 7”s. First set 
T = n and then a polynomial will be as follows: 


+ +H; + x2) 
+ x, + Xs) 


n+1 


XeXp-1 for odd n, k= 7 


+ Xe-1) for even n, k = 

+ [product terms of the degree equal to or more than the 
third, each of which includes at least one of the variables x; 
to xz, (for odd n, but to x, for even m) }. (59) 


That is, there exists a product term of the first degree, x,, a polynomial 
consisting of the second degree products that has some sort of regu- 
larity, and further for n 2 6, products of the degree equal to or more 
than the third. Now we shall determine when this polynomial con- 
sisting of the second degree products is invariant to permutation of the 
variables. In this polynomial the number of products, each of which 
includes x;, and the number of products, each of which includes x,_;, 


are jand n — j — 1, respectively, within a range of 1 < j < - D for 


odd n. For even n these are j and n — j — 1, respectively, within a 


range of 1 sj < — 1 and furthermore there are terms 


which include a single variable x,;2. Therefore this polynomial cannot 
be the same if we permute the variables, some of which have different 
numbers of products, but there are possibilities of obtaining the same 
polynomials after permutations of those variables that have the same 
numbers of products counted above. Such variables are x, and x,~1, 
n—1 
2 
nomial consisting of the second degree products in (59)) or x,_1 and 


obtained by setting j = for odd m (the last term in the poly- 


X~(= Xn/2), Obtained by setting 7 = ; — 1 for even n. Actually the 


following polynomials for » = 4, 5, and 6; x4 + X3(x2 + xi), 
+ x\X2X3, are kept invariant for permutations of x; and x2, x2 and x3, 
and x2 and x;, respectively. However for large m, if the third degree 


: 
hex 
: 
: 
: 
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terms in (59) are taken into account, this partial symmetry is broken, 
therefore giving different functions. That is, for odd m such that m 2 7, 
(that is, k — 2 > 1) x,.%,_2x, exists in the third degree terms including 
x, but the term x,_:%,—2x%, is not there, so we can discriminate by the 
third degree terms two functions which are made by permutation of x, 
and x,_;. Also for even ” such that nm = 10, (that is, k — 2 > 2) in 
the third degree terms including x, there exist two terms x,.%,-1%, and 
X-Xp_oX2 but the term x,_1X,-2X2 is not there. Then permutation of x; 
and x,_; does not give the same function. 

Consequently the number of functions given by permutation of n 
variables is 


8andn <6 


(60) 
= 7andn = 9. 


Next, decreasing 7 lower than n gives functions that are different 
from the above and one another. The number of the first-degree terms 
in (59) increases, that is, x7 + X74: + --: +x, for 7. With respect 
to variables x; to x7_1, terms corresponding to ones of higher than the 
first degree in (59) appear in the functions. Then permutation of n 
variables gives many different functions, the number of which is ob- 
tained by 


= 


for 7 =7and7 29 
and 


for T = 8and3 <T $6. 


Conversely increasing 7 starting at a value greater than m eliminates 
terms of the first degree. If we set J = n +s, a polynomial in the 
function, consisting of only the second-degree products, includes only 
variables from x, to x, for s such that 1 S$ s S$ m — 1 and then the form 
of this polynomial is identical with that in (59), that is, 


+ Xe-1), forevenn+s5,k (62) 


Also the function includes a generally complex polynomial consisting of 
terms of higher than the second degree. Though permutations of x, 
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and x,_; may give the same function, we can discriminate functions for 
n = 8 by taking into account terms of higher than the second degree as 
in the case of T = n. Then the number of different functions given by 
permutation of n variables is not less than 


n n! 
(63) 


for n = 8, where 1 Ss Sm—41. Further increase of 7 gradually 
raises the lowest degree of products in the functions and makes their 
functional forms more complex. 

Addition of three formulas (60, 61 and 63) over all possible 7 results 
in a number of the second type of functions given by permutation of 
variables. If m 2 9, division by 2 is not required for 7 29. A sum 


of (60) and (61), by use of (m — g)! (*) 


n! 
—., results in 
g! 


n! n! 
mits t (04) 
and a sum of (63) over s gives 
n! n!} 
nit + + (65) 
The sum of these is greater than 
1 
nt| 20 in 9)! (66) 


Calculations for cases of nm < 8 (some of which are shown in Table II) 
are possible but somewhat complex. The number of functions given 
by taking into account negation of variables of a number equal to or 
less than n can be obtained by multiplying (66) by 2" and adding them 
over integers not larger than n(2 8) but the number for gives a good 
evaluation. 

Last it may be shown as follows that the second type of functions 
differs from the first. In the second type of functions the polynomial 
consisting of only the second-degree products in each function has a 
partial symmetry, that is, it is invariant only to permutation of two 
specific variables (x, and x,-:). However, the first type of functions 
has one of the following cases: (1) a function does not have any second- 
degree product ; (2) or, if it has, the polynomial consisting of the second- 
degree product in a function is invariant to permutation of more 
than three variables when JT > 3. When JT = 2, we have xix; 
+ (x3 + --- +-~x,) which does not appear in any of the second type 
functions. 


= 
: 
d 
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However it should be noted here that even if repetition of some 
integers in a series of consecutive integers are repeated (like 433221), 
we still have functions with the same structure but with partial sym- 
metries. If we take account of this synthesis, evaluation of the bound 
could be further improved. Table II includes functions of this structure. 


Upper Bound on the Number of the Majority Decision Functions 


First we shall state the following important result known in the 
theory of linear inequalities. 


Lemma 1: Consider a set of consistent simultaneous linear inequali- 
ties (4)where w, ---, w, and T are regarded as variables. If the rank 
of its coefficient matrix is 7, this set of inequalities is satisfied by a set 
of solutions obtained by equating 7 inequalities in (4) (10). 


The following theorem is obtained with this Lemma. 


Theorem 16: All functions of variables that can be realized by a 
single majority decision element can be realized by properly selecting 
integers for coupling numbers of the input variables satisfying 


Proof: Assume that a given function is realized by a single ele- 
ment and then we have either the first type of an equation (4a) or the 
second type (4b), according to whether the function has the value 1 or 0 
for each of 2* sets of values of the variables. Then if this set of linear 
inequalities should be solved for (7+ 1) variables w,, ---,w, and T, 
Lemma 1 shows that if we choose and equate inequalities of a number 
n+ 1 out of the 2” inequalities, their solutions satisfy all other in- 
equalities. Writing these equalities as 


+0,—T 


(68) 


we have 


(69) 


Here A is a coefficient determinant of (68) where all elements in the 
first to the m** columns are either 1 or 0 and all elements in the (m + 1)** 
column are —1. A, is a determinant such that the 7** column of this 
determinant is replaced by the right sides of (68) consisting of 0’s and 
—1’s. Then divide all elements in the (7 + 1) column by 2. Next 
add these —} to other columns each of which consists of elements 1 and 


= 
2) 
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0, and subtract these —} from the column consisting of elements —1 
and 0. Thus all elements in the determinant 2A; are +} or —}. 
Application of the Hadamard’s evaluation of a determinant leads to 


(n+1)/2 
jad (2° (70) 


As A is an integer not smaller than 1, then w; = A,/A and 7 are rational 
numbers satisfying the inequalities. However if all w,’s and the 
threshold are multiplied by A, we have a logic element, with integral 
coupling numbers A-w;, which still represents the same function, 
because A,’s are integers. 

It may be interesting to compare the bound (67) with the largest 
coupling numbers required actually as follows. 


actual largest coupling 
number (16, 17) bound (67) 


NEDA 


1 
3 


NIA US 
VOM 


Theorem 17: An upper bound on the number of functions of vari- 
ables of a number equal to or less than m that can be represented by a 
single majority decision element is 


— 1) {4[” 41 (71) 


Proof: As in the proof of Theorem 16 it was shown that the maxi- 


n 
mum of A; is 2 [{ - , a number of combinations of integral 


coupling numbers ‘ae n variables is the n** power of this maximum. 
But as functions of less than n variables should be counted, w; = 0 for 
all 7’s should be included and then the number is at most 


(72) 


Since Theorem 6 showed that the number of possible values of T that 
give different functions is at most 2" — 1 excluding the two functions 
of 1 and 0 for each of the above combinations, the number of functions 
is at most the value given by (71), taking into account negation of 
variables. As this formula is counted with some redundancy, addition 
of two to this is not necessary in order to count the two functions of 


1 and 0. 
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In order to explain the extremely large memory capacity of the 
human brain as compared with the number of neurons, 10'*, D. G. 
Willis (11) proposed a hypothesis that variation of coupling numbers of 
inputs gives a single neuron the capability to represent various com- 
plex functions, providing it a large memory capacity. Quantitative 
evaluation for this may be obtained from the above two theorems. As it 
is said that a single neuron may have inputs as large as several hundred, 
if we assume n = 128, a lower bound on the number of functions repre- 
sented by a single neuron is 846 bits from (53), and its upper bound is 
about 4.1 X 10‘ bits from (71). For = 256, a lower bound is about 
1900 bits and an upper bound about 2 X 105 bits. 

From Theorem 16, we have the following Corollary. 


Corollary 4: All possible functions of » variables that can be realized 
by a single majority decision element can be realized with the integral 
values of total coupling number not larger than 


Proof: \f the threshold T is larger than > w;, we have only the 


constant function that takes the value 0 permanently. Hence, 
Ts tw, (74) 


should hold. As the relation (3) holds for positive integral coupling 
numbers of variable inputs and the maximum of the total coupling 
number also occurs with a constant input of 0, the total coupling number 
is 2T — 1 from (3). When w,’s are not integers, (74) and (67) hold 
even if A-w,’s are regarded as new coupling numbers. We therefore 
now have (73). 


5. CRITERION OF REALIZABILITY AND STRUCTURE DETERMINATION OF THE MAJORITY 
DECISION ELEMENT FOR A GIVEN FUNCTION BY LINEAR PROGRAMMING ‘ 

In order to realize a single majority decision element for a given 
function, we must determine the structure of the element so as to 
satisfy a set of the linear inequalities (4a and 4b) characterizing the 
function, but solutions cannot be uniquely obtained. For example if a 


‘ After preparation of the manuscript of this paper, the authors found two related papers; 
R. McNaughton, “Unate Truth Functions,’ Appl. Math. and Statistics Lab., Tech. Rep. 
No. 4, Stanford University, October, 1957 (appeared in Trans. IRE, Vol. EC-10, No. 1, 1960) ; 
R. C. Minnick, ‘Synthesis of Linear-Input Logic by the Simplex Method,”’ presented at the 
Sixth Annual Symposium on Computers and Data Processing of the Denver Research Institute, 
July, 1959 (ibid.). The authors believe that the present paper does not overlap with the first 
(except McNaughton’s theorem 3.2 in the Trans. JRE edition) and may include details different 
from the second, such as the reduction processes of the number of inequalities. 

Epitor’s Note: This manuscript was submitted to the JOURNAL office July 27, 1960. 
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set of solutions is discovered, when it is multiplied by an arbitrary 
number larger than 1, this gives another set of solutions. Then we 
shall add a restricting condition to minimize K = w,;+w:.+-::: 
+w, +w., the total coupling number of the element shown in (1). 
This restriction is actually important since it has economic and engi- 
neering merits. That is, the minimization of coupling amplitude 
among elements may result in a reduction of back couplings among 
elements and an increase of the number of outputs that can be obtained 
from a single element. These factors imply better stability of the 
circuits (15). When couplings among elements are realized by wiring 
as in the case of parametrons, this condition corresponds to minimizing 
the total sum of wiring turns. When elements are coupled with re- 
sistors as in a certain type of circuitry with the Esaki diodes, this corre- 
sponds to a maximization of resistor values. If it is desirable for the 
resistor values to be within a certain range, we must multiply a set of 
solutions by a number of proper magnitude or add inequalities repre- 
senting a range for the resistor values, to the inequalities (4a and 46). 

Adding the minimization condition to the set of the inequalities 
(4a and 46) in this way, determination of the structure of a majority 
decision element for a given function can be reduced to a typical prob- 
lem in linear programming that gives us at the same time a criterion of 
its realizability, as shown later. 

Since there are duality and inclusion relations between two majority 
decision functions, one with a constant input 0 and the other with a 
constant input 1, but with the same coupling numbers of variable and 
constant inputs as shown in the section on Duality of the Majority 
Decision Functions (p. 387), determination of the structure of the 
element can be simplified as follows: First check whether the given 
function can be realized by an element without a constant input, by 
finding whether the function is self-dual. Next, if this is not the case, 
compare the function with its dual and choose the smaller one (see 
p. 387). If it is realized by a single element, the element must have a 
constant input 0 and the other larger function can be realized by an 
element of the same structure except with a constant input 1. Conse- 
quently we shall henceforth only consider realizations of elements with 
no constant input or with a constant input 0. If there is no inclusion 
relationship between the given function and its dual, this function is 
not realizable, as seen from Corollary 2. Of course other various 
necessary conditions already discussed can be used to reject unrealizable 
functions. 

Now for simplicity, let us assume that all coupling numbers are 
integers and that K is odd. A more general case where some of them 
are non-integers will be discussed at the end of this section. 

As when all coupling numbers are integers the relation (3) between 
T and K can be used, the problem is simplified, excluding w.. 
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Problem 1: When all coupling numbers w,, and a threshold T to be 
found, are assumed as positive integers, determination of the structure 
of a majority decision element for a given function is reduced to a 
linear programming problem; the objective function, 7, is to be mini- 
mized under constraints of a set of the inequalities (4a and 4b). (The 
objective function is actually 2T — 1 but constants —1 and 2 can be 
ignored.) 

By the duality theorem in the theory of linear programming this 
problem can be converted into the following Problem 2 where the 
objective function is to be maximized. Solutions for this problem will 
be given by use of the simplex method and will simultaneously give 
solutions to the Problem 1. In this case there is no guarantee that the 
W;, We, -+-, W, and T obtained are integers, but these will always be ra- 
tional numbers since all coefficients in the set of inequalities are integers. 
When these happen to be non-integers, integral solutions can be ob- 
tained by the method developed recently by R. Gomory (12), but we 
shall not discuss this in this paper. 


Problem 2: Problem 1 is dual to the problem: under the constraints 


va; 0 (4 = 1,2, ---, m) 


j=l j=rt+l 


(75) 


j=l jert+l 


and 
(76) 


the objective function 


(77) 


j=rtl 


is to be maximized. 

The realization problem is thus converted into the regular type of 
linear programming problem, which can be solved by making tableaux 
of the simplex method. However since coefficients in condition (75) 
are only +1 and 0, it is generally a case of degeneracy and then in 
making tableaux the perturbation must be taken into account. By this 
perturbed version of the simplex method, if a given function is to be 
realized by a single element, w; (¢ = 1, 2, ---, m) and T can be deter- 
mined as optimum solutions and furthermore if it is not, that is, if 
Problem 1 does not have solutions, the non-existence of solutions is 
shown in the last tableau in Problem 2. 

The number of inequalities included in (4a and 40) is just 2", since 
there are » variables. Even for small values of » this number becomes 
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too large to be handled by electronic digital computers and solution by 
the simplex method may appear impossible. However, by the follow- 
ing processes we can cut it down greatly. 


Process A: Choose any one of the inequalities of (4a). If it is 
satisfied by variables w;, ---,w, and 7, any other inequality in (4a) 
such that any other terms are added to the chosen inequality, if such 
an inequality exists, is also satisfied. Then the latter inequality can be 
eliminated as a redundant one from (4a). Elimination of redundant 
inequalities in this way leaves inequalities each of which has the least 
number of terms. These will be called “the set of the minimal-term 
inequalities.” 

With respect to the inequalities of (40), if there exists an inequality 
each of whose terms appears in some other inequality that may itself 
have more terms than the original, then the original inequality can be 
eliminated as redundant from (4b). This elimination process of re- 
dundant inequalities leaves inequalities each of which has the largest 
number of terms. These will be called “the set of the maximal-term 
inequalities.” 

Now let us consider the Boolean lattice constructed on the partial 


ordering defined by 
(78) 


for two sets, X and X’, each of which consists of numbers from 1 to n. 
As easily seen, the maximal-term inequalities correspond to the gen- 
erators of the principal ideals in this Boolean lattice and the minimal- 
term inequalities correspond to the generators of the principal dual ideals. 
Each element of the lattice must belong to one of these ideals or dual 
ideals. 

The process of extracting the minimal-term inequalities from those 
of (4a) is equivalent simply to an application of the deletion law de- 
scribed in the paragraph after Definition 2, that any term including 
some other term as its factor in the function expression can be deleted. 
Then we can write down the minimal-term inequalities directly from 
the given function as each term in the polynomial expression for the 
function, that is, each prime implicant has one-to-one correspondence 
with each of the minimal-term inequalities. Now there is a question 
of whether we can write down the maximal-term inequalities from this 
polynomial. The procedure is as follows. 


Process B: Make the polynomial dual to the given polynomial and 
then eliminate all redundant terms by application of the deletion law. 
Next, for each term included in this dual polynomial make a product 
consisting of only variables that are not included in the term, and now 
all such products specify the maximal-term inequalities, each product 
corresponding to each inequality. 


: 
: 
: 
4 
: 
3 
is 
i 
; 
; 


406 S. Muroca, I. Topa anv S. Takasu (J. F. 1 


This is easily seen from the fact that negation of the original func- 
tion specifies the maximal-term inequalities and it is equivalently ex- 
pressed by its dual function with negated inputs. 

For example, when the function + x1%3 + + is given, 
the minimal-term inequalities corresponding to each term will be as 
follows, 


W, + We —T20 
w W; -—-T20 
1 + ws = (79) 
Wi + th T = 0 


Figure 1 shows the Boolean lattice of the minterms where the lattice 
elements are indicated by suffixes of the variables. The inequalities in 


Lattice for four variables. 


Fic. 


(79) correspond to the lattice elements with circles. Next, as its dual 
function is + + according to the process above we 
have a sum of three terms x3X4, X2X, and x; corresponding to the respec- 
tive terms of the function. The lattice elements with the triangles in 
Fig. 1 indicate these terms and, together with the lattice elements with 
the circles, bisect the whole lattice, where ¢ indicates the lower bound 
of the lattice. The maximal-term inequalities are obtained from these 
terms as follows: 


—w,—wuw,+T7 21 
— Ws (80) 
— WwW, +721. 


123 124 134 234 
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These processes reduce 2‘ = 16 inequalities to only seven inequalities in 
(79) and (80). 

The number of these inequalities may be, in general, reduced further 
by the following process of comparing coupling numbers of variables 
indirectly, though of course we shall get exactly identical solutions 
without this process, if solutions exist. 


Process C: Choose any two terms A and B in the given original func- 
tion and compare them in the following way. Let f4,» be the function 
obtained by assigning values of 1 to the variables included only in the 
term A of the original function and values of 0 to the variables included 
only in the other term B. This is a function of the residual variables. 
Next, let fz,4 be the function obtained by assigning 0 to the former 
variables and 1 to the latter. If f4;2 D fx), eliminate an inequality 
corresponding to term A from the minimal-term inequalities and, if 
fare Chea, eliminate an inequality corresponding to B. If neither 
holds, do not eliminate either. Repeat this process over the remaining 
terms. Finally, when any inequality cannot be further eliminated for 
selection of any two terms, make inequalities corresponding to these 
terms remaining. In many cases, this process eliminates further in- 
equalities from the minimal-term inequalities. 

The reason why elimination of inequalities is possible by this process 
is explained by the following considerations. A set of feasible solutions 
of Problem 1 constitutes a bounded polyhedron and the optimum 
solutions correspond to one of its vertices. Then the optimum solu- 
tions are obtained by solving a set of simultaneous linear equalities 
obtained by removing the greater-than condition from the inequalities 
chosen properly out of the set of the inequalities (4a and 46) (13, 14). 
However, fa;z > fx;4 means that a sum of coupling numbers of vari- 
ables included only in A is larger than that included only in B. Then 
as it means that the left side of the inequality corresponding to A is 
larger than that corresponding to B, if only an equality “‘=’’ is to hold in 
the latter inequality ‘‘=,’’ only an inequality ‘‘>”’ but not an equality 
“‘="’ should hold in the former inequality. Consequently the former 
inequality corresponding to A is unnecessary, leaving only the equated 
constraint corresponding to B in the set of the simultaneous linear 
equations. However it is to be noted here that since we cannot in 
general eliminate all of the redundant inequalities, we do not know which 
of the remaining inequalities should be equalities and consequently we 
must solve a set of inequalities, as 2 and not =. But when the func- 
tion is a function of truly m variables, that is, when there are just 
n + 1 constraints, all of them can be equated. 

We shall explain the above process on the function used earlier: 
+ + + First compare and x2x;. A variable 
included only in the former is x; and one included only in the latter is 
x3. A function obtained by setting x; = 1 and x; = 0 is x, + x4 and 
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one obtained by setting x; = 0 and x; = lis x. As x2 + %4 D X2, an 
inequality corresponding to x,;x, can be eliminated. Next compare 
x\X3 and xeox;. A function given by setting x, = 1 and x, = 0 is 
x3 +x, and one given by x; = 0 and x, = 1 is x;. As x; 
holds, an inequality corresponding to x,;x; can be eliminated. As to 
the residual two terms, xx; and x:x;, such comparison of functions is 
not possible and then only inequalities corresponding to these two 
remain. 


Process D: Process C should be again applied, comparing products 
obtained in Process B, except that inequalities corresponding to term B 
instead of A should be eliminated from the maximal-term inequalities, 
if faye D feja. The reason for this is the same as in the previous case, 
except that coefficients of w,’s are negative. 

For the above example, Process D should be applied on the original 
function + + + By Process B, we had x3x4, 
and x; For example, compare xx, and x3x4 A function given by 
setting x. = 1 and x; = 0, in the original function, is x, and a function 
given by setting x, = 0 and x; = 1 is x, identical to the first function. 
Also, for any two terms, we have the identical function. Then no 
inequality should be eliminated. 


‘Process E: The above processes covered Problem 1. We must 
consider the dual Problem 2 with the eliminated set of inequalities after 
these processes. Then the number of inequalities in (75) is m + 1 but 
the number of variables v,’s is in general greatly reduced from 2". 


A Numerical Example: Applications of the above Processes A—E 
on the inequalities (4a and 4) for the function x\x. + x1x3 + xox; 
+ x:x, reduce 16 inequalities to the following five : 


This is converted by Process E into the problem of finding the optimum 
values for v7, ---, ¥s such that under constraints 
Ve — Us 5 
Ve — U3 — 
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TABLE I.—Simplex Tableaux. 
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and 
(83) 


the objective function 


V3 +45 (84) 


should be maximized. 

Now proceed with the solution of this problem using the simplex 
tableaux, taking into account the perturbation, as shown in Table I. 
Finally, if the problem has solutions, the structure of the majority 
decision element to realize the given function will be obtained. In this 
table the row “‘u;’’ shows coefficients of the terms in the objective func- 
tion (84). Andz; = >-u;b;; where is an entry in the column and 


the j** row. Ae, and are slack variables and the column “‘c,” 
shows the right-side constants of the constraints. 

As shown in the last row of the optimum tableaux, the optimum 
solutions that minimize the total coupling number as 27 — 1 = 7, are 
W, = 3, Ww. = Ws = 2, ws = 1 and T = 4, and this structure certainly 
represents a given function. 


ProcessgF: The tableaux are made as follows. First choose any 
entry, 2 — u,, with a negative sign and with the largest absolute value 
in the row labeled as ‘‘z; — u;,’’ and then find },:.’s with a positive sign 
in the column to which the above z — u, belongs. Since generally 

1 
there are several b,,’s, calculate (Ca + b.;€’) for each of such b,,’s, 
where c, is the entry in the column labeled as ‘‘c,’’ and where the sum- 
mation > is over the entries in the columns of the slack variables 


(j = 1,2, ---, 5 in Table I) and where « is a sufficiently small positive 
number. We can find only one d.+. such that the above formula has 
the minimum value. Next, replace all entries }.+; in the a*** row by 
as; = except that the entry in the column labelled “variable,” 
and in the a*t* row, is replaced by the symbol in the k** column and in 
the row labelled ‘‘variable’’ and that the entry in the column labelled 
“u,;’ and in the a** row is replaced by a number in the row labelled 
“u;’’ and in the k** column. Then, replace all the other entries in the 
k* column by 0 and the other 0;;’s in the tableau by );; — b’a«;-dix. 
Finally, calculate z; — u;. Thus, the next tableau was obtained. 

For the next tableaux repeat the above process which is one usually 
used in the simplex method. 

In the optimum tableau we shall finally have the following three 
cases. Here it is important to note the following: If we apply the sim- 
plex method on a set of inequalities reduced by Processes C and D, 
there is no guarantee for solutions obtained in the following cases (i) 
and (ii) to be solutions for the original Problem 1, too. Therefore we 


i 
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must check whether the solutions obtained satisfy a set of inequalities 
obtained by applying only Processes A and B, or alternatively, whether 
the coupling weights and threshold obtained construct a given function. 
If there is no solution as in the following case (iii), the original Problem 1 
also has no solution. However, when the simplex method is used after 
applying only Processes A and B, no such checking procedure is required. 


(i) Consider the columns of variables (the 6 to 10* columns in 
Table I). If z; — u; > 0 holds for all columns except ones (adopted 
as the bases) where all entries except a single one are replaced by 0 in 
the above process, we have a single set of the optimum solutions for 
Problem 1 given by (z; — u;)’s in the columns of the slack variables. 

(ii) If z; — u; = 0 holds for some columns (in columns other than 
the bases), we have other sets of the optimum solutions than the one 
in the optimum tableau. 

(iii) If we cannot find any b,, with a positive value in Process E, in 
other words, if only zeros or negative values are found in the column, 
the solutions for Problem 2 are called ‘‘infinite’’ because the objective 
function can be made arbitrarily large. In the original Problem 1 this 
means that there exist no solutions. When a given function cannot be 
represented by a single majority decision element, the last tableau 
always takes this form. 


In the last tableau, if some slack variables remain in the bases, 
coupling numbers corresponding to them become 0. This is a case 
where a given function can be represented with fewer variables than 
those apparently required. It is interesting to note that if the coupling 
numbers and the threshold are not 0, the inequalities corresponding to 
these in (75) should hold with equalities “‘=” only (14). 

So far no guarantee has been obtained that solutions w;, ---, w, and 
T obtained by the above simplex method will always be integers. 
However, as mentioned later, we had observed that the optimum solu- 
tions w,’s and T are integers for any functions with variables of a small 
number n, though non-integers appeared in intermediate tableaux in 
some cases. 

In the general case where there is a possibility that such non- 
integers may be included among coupling numbers of the inputs, it 
may not be reasonable to solve Problem 1 directly, where the objective 
function T was set up under the assumption that all coupling numbers 
are integers. Then we must solve the problem where the objective 
function K = w,;+--- +w,+w. should be minimized under the 
constraints (4a and 4b). However, in this case, T is included only in 
the constraints while w- is included only in the objective function. Then 
a relation among 7 and wi, w2, --:, w, and w., should be added to the 
constraints. To avoid this, let us start from the following expression. 
An equivalent of (2a and 2b) in Definition 2 of the majority decision 
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elements, where w, is used instead of T is: the output value of the 
majority decision element is 


14 + Ox. ~"Dw, >'0 


i=1 


O, if — 1) + (2x, — 1m, < 0, 


i=1 


where x;;, «++, X,; are components of the j** input vector (7 = 1, ---, 2") 
and x, is the value of the constant input. 


Problem 3: Determination of the structure of a majority decision 
element for a given function is reduced to the linear programming 
problem: the objective function K = w;+-:-: +w, +, is to be 
minimized under constraints of a set of inequalities (85) where equali- 
ties ‘‘=”" are included in all inequalities after knowing the minimum of 
the left sides of the first inequalities and the maximum of those of the 
second. 

Here it is important to note that if the inequalities in (85) do not 
include equalities ‘‘=,’’ then arbitrarily small w,’s and w, can be solu- 
tions. This is an essential point in this minimization problem. Then 
the minimum .¢, of the left sides of the first inequalities and the maxi- 
mum €¢2 of those of the second inequalities should be known or pre- 
determined. These give a margin in the performance of the element. 
Choose the smaller one out of ¢€; and |¢:| and write it as ¢, that is, 
e = Min (e:, |es|). All inequalities in (85) can be made to include 
equalities “‘=,”’ by normalizing the inequalities by dividing these by e. 
Therefore, Problem 3 may be used instead of Problem 1 and the simplex 
method can be applied again. Exactly the same processes as those 
used on Problem 1 can be applied here, as easily seen. This form of the 
problem may be more basic since we are working directly on w, without 
T, but may be somewhat complex. : 

However, it is important in this case to see that if we introduce a 
parameter T characterizing the majority decision element given by 
T = 3(K +1) after this normalization, the constraints in Problem 3 
coincide with those in (4a and 46) and the objective function also will 
be T. Consequently, Problem 3 coincides formally with Problem 1 
and then even in the general case where non-integral coupling numbers 
may be included, it may be sufficient for us to solve Problem 1 formally. 

Here it should be noted that in Problem 1 the relation between T 
and w., (3) is not included. The minimization of the threshold T only, 
as the objective function, without taking into account relation (3), 
implies intuitively minimization of sums of coupling numbers which 
do not reach the threshold. 
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6. MAJORITY DECISION FUNCTIONS OF UP TO SIX VARIABLES 


In order to check various properties of the majority decision func- 
tions described so far, the Parametron Computer MUSASINO-1 (15) 
was used to find all majority decision functions of up to six variables 
and the structures of elements to realize these, based on Problems 1 and 2. 


Majority Decision Functions of a Few Variables and Their Number 


All majority decision functions of up to six variables are classified 
by permutations of variables, and negations of the variables. A 
representative function in each class such that w,; 2 w, = w; = wy 
> ws = we = O are obtained (16, 17). The number of these functions 
is too great to be shown here. 

These functions were checked both by the simplex method and the 
combination method, where the functions are discovered by giving all 
possible integer values to the inputs. 

Table II lists the numbers of functions enumerated from several 
different viewpoints (18). 

(i) The number of general Boolean functions of up to m variables, 
22", is shown for reference in the second column of Table II. 

(ii) The number of the positive functions of up to n variables, 
those that can be realized without negation of the variables. These 
may generally not be realizable by a single majority decision element. 
The numbers were quoted from Birkhoff’s book (19). The majority 
decision functions constitute only a subset of these. 

(iii) The number of the majority decision functions of exactly n 
variables without negations of the variables. This is divided into two 
cases. The first case shows the number of the majority decision func- 
tions, representatives of equivalent classes, each of which consists of 
functions identical by permutations and negations of the variables. 
For example, x1%2 + x3, x1%; + x2 and #::%: + x; belong to the same 
class represented by x; + x2x;. The second case shows the number of 
the functions in the first case, taking into account only permutations of 
variables. The functions in the second case are not enumerable unless 
the forms of the functions representing the above classes are known, 
because there are partially or totally symmetric functions. For ex- 
ample, for m = 2, we have two, counting x; + x, and x;%>. 

(iv) The number of the majority decision functions of up to m vari- 
ables, without negations of variables. For example, for = 2 we have 
six, counting 0, 1, x1, %2, x1 + x2 and Xx1%2. 

(v) The number N(m) of the majority decision functions of up to 
variables, taking into account permutations and negations of variables. 
If this is compared with 22", we can see how small a part of 2%” this 
occupies. 

(vi) Theoretical upper and lower bounds on N(m) shown in Section 4. 
The lower bound shows the number of functions constructed similarly to 
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the methods in Theorems, 8, 13, 14, and 15. (Cf. the paragraph just 
preceding Upper Bound on the Number of Majority Decision Functions, 
p. 400.) As the bounds in Section 4 were aimed for cases of large n, 
these approximations are not good here. 

The ratio of the number of majority functions to that of the general 
functions or to that of the positive functions approaches zero as 
increases. 
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x) +X2X3 
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Lattices of the Majority Dectston Functions 


Though there may be various ways to classify the majority decision 
functions, a classification based on the ordering among the input 
variables and that among the input vectors in Section 3 will be dis- 


cussed here. 


(1) The lattice based on ordering among the variables. 


If an ordering x; > x; between any two variables x; and x; defined 
for a function f;(x1, - X,) holds also for another function f2(x1, Xn), 
we can define an ordering between x; and x; for the functions f, + f2 and 
fife where x; > x; still holds. Therefore, the whole set of functions 
that defines ordering among the input variables, for example, 
2X2 2 Xn, Clearly constitutes a distributive lattice with 
respect to the operations ‘‘+’’ and “-.” This lattice will be denoted 
as L,”. 

Though functions constituting the above lattice are not always the 
majority decision functions, it was ascertained that such functions are 
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confined to the majority functions only when the number of variables 
does not exceed five. Then the whole set of majority decision functions 
such that w; = 2 = w,, coincides with L,“ for n s 5. 

The Hasse diagrams of the L,’s for m < 4 are shown in Figs. 2 
and 3. As shown in Theorem 5, when we increase the threshold of 
an element with fixed coupling numbers of variable inputs, the functions 
get smaller, as fo Df: D --- Dfw, where f; is a function with the 
threshold z. A series of such functions fo, fi, -- -,f is a subset in a chain 
of the L,“. Also, from Corollary 2, this subset proves to be sym- 


metrical with respect to the straight line S in Fig. 2. 
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1 
| 
+xX2+xX3 


xj) +x2+x3 


x) +xX2 
| 
+X2x3 


x) +X2x3 


(x2+x3) 


x} (x2+x3) 


%1%2 *1%2 
| 


*1%2%3 


0 0 

f f f f f f 
(111) > (110) > (101) > (100) (111) > (110) > (101) > (011) 

f f f f ee f f 
> (011) > (010) > (001) > (000) > (100) > (010) > (001) > (000) 


Fic. 4. Examples of L;®). 


(2) The lattice based on orderings among the input vectors. 


Similar to the above, if any two input vectors Z; and =; are com- 


parable simultaneously for two functions f; and fs, as Z; > 2; and 


fitfs 
=; > Z;, then also Z; > Z; and =; > Z;. Therefore, the whole set of 
the functions which arrange all the input vectors {Z,} in a certain order 
proves to constitute a distributive lattice with respect to the operations 
and 

We can construct various lattices according to what ordinal relations 
are given among {Z;}. However, we cannot establish the relations on 
{Z,} in a completely arbitrary way so that these can be defined by 
some functions. 


Now let L,® be a special lattice such that orderings 2; > 22 
f 
> --++ > Zs are established and x; 2 x. 2 -:- 2x, hold. L,® isa 


sublattice of Z,“. It was ascertained in a range of m <5 that all 
functions included in L,® are majority decision functions and some 
chains of 

Figure 4 shows examples of L;‘*. These exhaust all chains of 
L; but in the cases of m = 4 and 5, not all chains of L,“ are included 


in L,®. 
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READER OPINION 
SPIRAL STELLAR NEBULAE AND COSMIC GAS JETS 


BY 


C. E. R. BRUCE! 


Since I sent you my paper on the above subject (1)? | have had my 
attention called to O. C. Wilson’s paper (2) on ‘The Nuclear and 
Nebular Spectra of the Planetary Nebula NGC 2392.” It will be 
recalled that this is the nebula which Campbell and Moore (3) took as 
the type nebula of the group in which the gas movements proved quite 
inexplicable on any of the many hypotheses they tried. I showed in 
my Fig. 4 that these movements were at once resolvable into the two 
diametrically opposed gas jets, initiated by the two electrical discharges 
which had occurred at some time in the star’s atmosphere. The tem- 
perature, and hence the gas velocity will vary along each of the two jets 
in a more or less similar manner, so that any particular level of excita- 
tion will be found at corresponding positions in each of the two jets. 
Thus if a slitless spectrogram of the nebula is obtained it is to be ex- 


Photograph by courtesy of Mount Wilson Observatory 


pected that there will be two rings, one displaced to the violet, and one 
to the red of the normal position of the line. This will be seen in the 
accompanying figure to be confirmed by the slitless image of the line 
43868 [NellII], taken from Wilson’s paper [his Fig. 2b). The brighter 
of the two rings is the red component of the line, indicating that the 


1 Electrical Research Association Laboratory, Leatherhead, Surrey, England. 
? The boldface numbers in parentheses refer to the references appended to this note. 
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much greater intrinsic brightness of the further discharge more than 
counterbalanced the effect of the increased absorption of the light in 
the star’s atmosphere, as compared with that experienced by the light 
from the nearer of the two jets. Curiously Wilson refers only to the 
inner “ring,”’ and ignores the existence ofthe fainter neighbor displaced 
towards the violet. The outer ring from*the nearer jet is evidently too 
faint to be photographed, and only that associated with the ‘‘red”’ 
inner ring is seen, itself very much fainter than the corresponding inner 
ring. The presence of the level of excitation corresponding to \3868 
[NelIII] will depend on both the temperature and the density of the 
gas. The outer ring is apparently emitted when the jet has almost 
reached the limits of the star’s atmosphere, where presumably the 
density and the temperature have both decreased, but in such a manner 
that their combination yields conditions favorable to the existence 
of the same level of excitation as appears lower down in the star’s 
atmosphere. 

However, despite the absence of the violet component of the outer 
ring, the appearance of the two inner rings confirms the existence of 
the two jets depicted in Fig. 4 of my paper. 
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A NOTE ON THE DUCTILITY OF BERYLLIUM SINGLE CRYSTALS 
ORIENTED FOR BASAL SLIP AND TESTED IN TENSION* 


BY 
M. HERMAN AND G. E. SPANGLER 


This note has been prepared to present experimental evidence which 
suggests that impurities in beryllium strongly influence the amount of 
basal plane glide obtained during tensile testing. 

Several single crystals of beryllium were prepared by zone melting 
using the floating zone method. The apparatus used was similar to 
that previously described for the zone refining of titanium (1)! with, 
however, the addition of a gas re-circulating pump to limit the vapor 
deposit in the vicinity of the molten zone. Crystals obtained from 
this process were machined to conical shouldered tensile specimens by 
electro-chemical methods (2) or a combination of spark-discharge and 
electro-chemical machining. With the exception of one bend-test 
specimen (see Fig. 1) the specimens were tensile tested in an Instron 


Fic. 1. A tensile specimen (prepared from 8 “‘pass’’ material) which was subjected to 
bending. Orientation, xo approximately 25°; bend axis [1010]. Elongation on outer fibers 
approximately 50 per cent. (Magnification, 1.2 X.) 


Testing Machine in which the load-extension curves were automatically 
recorded. 

In Table I the history of the specimens and the tensile testing results 
are presented. 

These results, in contrast to previously published work (3) on the 
tensile deformation of single crystals of beryllium, clearly indicate that 
the extent of basal glide occurring during tensile testing is not necessarily 
limited to very small values. Based upon the apparent relationship 
between the extent of zone melting (number of passes) and the values 
obtained for the resolved shear stress at yielding, it is suggested that the 
observed increase in ductility is the result of a decrease in the impurity 
concentration in the beryllium. 

These data are being examined in terms of current theory of the 


* This work was supported by the U. S. Navy, Bureau of Naval Weapons. 
1 The boldface numbers in parentheses refer to the references appended to this note. 
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TABLE I.—History and Tensile Testing Results of Beryllium Single Crystals. 


Starting 
Material 


Pechiney (3) 


Berylco 
Vacuum Cast 
Ingot 


Berylco 
Vacuum Cast 
Ingot 


Pechiney 


Number of Zone 
Melting Passes 


Vacuum Cast 


Heat Treatment 
rior to 
Testing 


(information 
not available) 


Orientation* 
xo=45° 

Xo (infor- 
mation not 
available) 


xo=47° 800°C.-1 hr. 


Ao 51° 


xo = 20° 900°C.-1 hr. 


Xo = 25° 


xo =48° 
Ao = 55° 


* xo = angle between slip plane and tension axis 
do = angle between slip direction and tension axis 


Resolved Strain at 
Shear Stress Fracture, 
at Yielding, % Glide 

psi. Strain 


2000 + 500 


Operating 
Glide 
System 


Basal Plane 


16 Basal Plane 


64 Prism Plane 


156 ~=Basal Plane 


cleavage fracture of single crystals while further experiments are being 
conducted to examine the details of the proposed impurity effects. 
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NOTE ON THE EFFECTIVE CONSTANTS OF 
COMPOSITE MATERIALS * 


BY 


Z. HASHIN! AND S. SHTRIKMAN ? 


The present work deals with the determination of effective material 
constants of two-phase isotropic materials, assuming that the macro- 
scopic physical behavior of each of the constituents is analytically 
known. The specific problem treated is that of the magnetic perme- 
ability of two-phase materials. For reasons of mathematical analogy 
the results obtained also apply to dielectric, electric conduction, heat 
conduction and diffusion properties of such materials. 

The problem is one of long standing (1)* and many different results 
have been obtained. It has been shown by Brown (2) that the effec- 
tive magnetic permeability is not completely defined by fractional 
volumes of the constituent phases and that statistical aspects of phase 
dispersion influence the results. It is the purpose of this work to find 
the most restrictive upper and lower bounds on the effective constants 
that can be obtained in terms of phase constants and fractional volumes, 
assuming isotropy of the composite material. 

These bounds are constructed by using the following variational 
principles (3). Consider a body of volume V and surface SS of material 
of magnetic permeability uo, which for simplicity is taken to be constant 
in space. Let a magnetic potential yo = Yo(.S) be prescribed on the 
boundary. This defines uniquely the field intensity in the body which 
is denoted by Ho. 

Consider now the same body where part or whole of it has been 
changed to material with permeability 4, which may vary in space. 
Define the generalized magnetization by 


T= 5 - poll (1) 


where B is a magnetic induction and H a magnetic field intensity. 


* The authors gratefully acknowledge support by the Office of Naval Research and by the 
Air Force Office of Scientific Research of the Air Research and Development Command. 

1 University of Pennsylvania, Philadelphia, Penna. On leave of absence from the Tech- 
nion-Israel Institute of Technology, Haifa, Israel. 

The Franklin Institute Laboratories for Research and Development, Philadelphia, 


Penna. On leave of absence from the Weizmann Institute of Science, Rehovot, Israel. 


5 The boldface numbers in parentheses refer to the references appended to this note. 
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Define 
H =H Hy. (2) 


Then the volume integral 


27H) dV (3) 


( 
Ur Sr pollo ito 


is stationary for 
T = (u — wo)H (4) 


where #7’ and 7 satisfy the subsidiary condition 
uo div H’ + div T = 0. (5) 
The intensity H’ must satisfy the boundary condition that y’, defined by 
H’ = — grad y’, (6) 


vanishes on the boundary S. 

The expression Ur is an absolute maximum when yp > wo and an 
absolute minimum when yw < wo. The stationary value of U7 is equal 
to the magnetostatic energy stored in the changed body. 

A similar variational theorem has been established for a body with 
prescribed normal component of boundary induction. 

Let a two-phase isotropic and quasi-homogeneous material of volume 
V consist of two phases of volumes V; and V2. The fractional volumes 
of these are denoted by v; and v2 where 


V1 + ve = 1. (7) 
Let yw: and ue be the permeabilities of the phases where 
Ho > (8) 


and 7, and 7; the generalized magnetizations in them. 

The effective permeability u* is defined as follows: Let a potential 
Wo, derived from a uniform H, be applied to the boundary of the com- 
posite body. The magnetostatic energy U* of the body is then taken 
to be given by 


1 
U* = (9) 


: 
q 

q 3 

j 
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With reference to Eqs. 1-5 it is now assumed that 


Mo = M1 
B= He (10) 


i= 0 
T. = Ty constant and parallel to Hp. 


(11) 


Introducing (10) and (11) into (3) and using (5), the following result 
is obtained 


T0102 
Pi 


Ur= ; + 2T Ve (12) 


The third term in (12) which corresponds to the third term in (3) 
requires some clarification. It has already been given by Néel (4) and 
can be derived by solving (5) subject to (6) by Fourier methods (5), or 
by making use (6) of the Brown-Morrish (7) theorem. 

Introducing (9) into (12) a lower bound on y* is obtained. The 
highest lower bound is found when Uy, in (12) is maximized. The 
result is 


* v2 


Me — 


To find an upper bound it is assumed that 


(14) 


T, = T» constant and parallel to Ho 15 
= 0. 


By the same procedure it follows that the lowest bound is given by 


* 
+ i (16) 


Mi — 


It should be noted that these bounds are independent of the geometry of 
structure of the composite material, the only requirement being isotropy. 

It can furthermore be shown that these bounds are the most restric- 
tive that can be obtained in terms of the phase constants and fractional 
volumes. This is done by proving that for certain phase configurations, 
ui* and ps* given in (13) and (16) are exact expressions for u*. For this 
purpose consider composite spheres each of which consists of an inner 
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sphere of one phase surrounded by a concentric spherical shell of the 
second phase, the volume ratio of the two phases being fixed in all 
spheres. Let a body be completely filled with composite spheres (this 
is possible if no size restriction is imposed). In this case, when the 
inner spheres have permeability wu. and the shells permeability yu; 


* (17) 


If the materials are interchanged, then yw.* = u*. This follows by a 


method used by Hashin (8). 
The theoretical results given in the present work are verified by 
experimental investigations of DeLoor (9) for the analogous dielectric 


problem. 
Similar results have been derived by the authors for two-phase 


elastic materials, formulating and using variational theorems analogous 


to the ones given here (10). 
The authors are indebted to Dr. W. F. Brown, Jr. for valuable ideas 


generously communicated. 
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COMMITTEE ON SCIENCE AND THE ARTS 
(Abstract of Proceedings of Stated Meeting held Wednesday, April 12, 1961.) 


HALL OF THE COMMITTEE, 
PHILADELPHIA, APRIL 12, 1961. 


Mr. Francis G. TATNALL in the Chair. 
The following reports were presented for final action: 
No. 3362: Work of Reinhold Riidenberg. 


This report recommended the award of an Elliott Cresson Medal to Reinhold Riidenberg, 
of Belmont, Massachusetts, ‘“‘For his many contributions to the electric power industry for 
over fifty years as an engineer, inventor, author and teacher, and particularly for his scholarly 
textbook on transient performance of electric power systems by which present and future 
engineers may attain a better understanding of the problems of their profession.” 


No. 3384: Negative Impedance Booster. 


This report recommended the award of an Edward Longstreth Medal to Josiah L. Merrill, 
Jr., of Murray Hill, New Jersey, ‘For his analysis of problems associated with losses and dis- 
tortion in telephone transmission systems; for his synthesis of the negative-impedance booster ; 
and particularly for his successful development of a negative-impedance booster system that 
has proved to be highly effective and economical.” 


D. S. FAHRNEY, 
Secretary to Committee 
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Sensory DEPRIVATION, edited by Philip Solo- 
mon, Philip E. Kubzansky, P. Herbert 
Leiderman, Jack H. Mendelson, Richard 
Trumbull and Donald Wexler. 262 pages, 
diagrams, 54 X 8} in. Cambridge, Har- 
vard University Press, 1961. Price, $5.75. 


Man is always searching for answers to the 
mysteries about his kind. And especially 
about his central function, which puts him 
at the top of the present phylogenetic tree. 
When the time comes that we can roadmap 
the anatomy of the brain, spinal cord, and 
peripheral nerves, down to the last neuron, 
maybe we can reach some certainties. With 
our present techniques and human limitations, 
this is too great a task. 

The physiology of single nerve function and 
transmission has been examined and under- 
stood much more thoroughly than the anat- 
omy or physiology of the central nervous 
system as a unit. In 1957, the mechanism 
of transmission of a single impulse from the 
endoplasmic reticulum to the cell membrane 
was demonstrated by the British. Only in 
recent months have the neurohormones and 
neuroreceptors involved in transmission be- 
tween cell membranes been elucidated. Our 
understanding of central function lags far 
behind this. 

The problem of collecting objective data 
from a system so beset with undefinable and 
unknown variables is truly monumental. 

In the last thirty five years techniques 
have been devised for obtaining data along 
many parameters of central function. Out 
of these developed several useful hypotheses. 
But the limiting problem is always the same: 
in the infinity of variables, there is not one, 
or preferably several hundred thousand, 
which can be adequately controlled and in- 
terpreted. Were this possible, we would have 
a chance to develop a starting point. Of 
course, there are anatomical and ethical 
prohibitions to objective measurement, too. 

The real problem in initiating a systematic 
study of central function is breaking into the 
circle of variables with one constant. Al- 


though the literature contains many repro- 
ducible observations, this goal is probably 
some years in the future. Studies of sensory 
deprivation may shorten this interval. 

Stemming from the investigations of “‘brain- 
washing”’ during the Korean War, and carried 
forward with mushrooming impetus in an- 
ticipation of manned space travel, the concept 
of depriving a human of the common stimuli 
to which he is constantly exposed (light, taste, 
smell, sound, light touch, deep pressure, heat, 
cold, pain, the kinesthetic senses, etc.), and 
studying the responses to this new stimulus 
deficit, has provided an entire new area 
for behavioral investigation. The field is 
only several years old, and so far, not much 
common ground has appeared among 
investigators. 

The purpose of this symposium in 1958 was 
to bring some unity and standardization to 
the area. As well as offering some excellent 
criteria toward this end, a cross section of 
revealing behavioral experiments are unfolded, 
adding to our understanding of why you and 
I do what we do. 

From the data many useful conclusions 
are evident. New versions of the stimulus- 
response hypothesis are offered, and the new 
concept that the reticular substance of the 
brain stem has cognitive and integrative 
functions is suggested. 

Though this is all valuable in itself, the 
measure of greatness of this investigative tool 
is that it may provide the break in the circle 
—the first controlled variable—not by adding 
measurable stimuli to the system, but by 
removing them from it—so that we may begin 
the systematic study of central function. 

This symposium is well balanced enough 
to be considered a standard textbook of 
psychology at the graduate level. The ma- 
terial is of interest to psychiatrists, and the 
greater meaning of it will be challenging to 
all biological philosophers. 

CHARLES W. ANDERSON 
Jefferson Medical College 
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QuANTUM THEORY OF ATOMIC STRUCTURE, 
Vor. II, by John C. Slater. 439 pages, 
6X9 in. New York, McGraw-Hill Book 
Co., Inc., 1960. Price, $13.00. 


While this second volume of Quantum 
Theory of Atomic Structure (Chapters 17 to 
26) is a natural sequence of the first volume, 
the level of treatment is considerably higher 
and the subject matter much more specialized. 
The Hartree self-consistent field method was 
dealt with in the first volume. The second 
volume begins with the Hartree-Fock method. 
The treatment of this method is thorough 
and beautiful, and many recent contributions 
of Slater and others are included, in particular 
an illuminating discussion of the “exchange 
potential.”” After a discussion of correlation 
problems, the application of group theory to 
atomic structure is exposed clearly but rather 
briefly (30 pages). The quantum mechanical 
tools thus sharpened are used to study more 
rigorously the more complex non-relativistic 
aspects of atomic spectra, namely, Multiplet 
theory and Racah’s coefficients. Dirac’s rela- 
tivistic theory of the electron follows, and 
its physical implications are fully discussed. 


The magnetic behavior of complex atoms, and 
transition probabilities of radiation are anal- 
yzed. The last chapter is devoted to hyper- 
fine structure (magnetic dipole, and electric 
quadrupole H.F. structure of one electron 
atom). 

The book is carefully and pedagogically 


written. It is an up-to-date synthesis of the 
subject matter considered, since so many of 
the contributions of the last twenty years have 
been interpreted for the first time in a single 
coherent picture. Some of the appendices 
and tables will prove to be considerably help- 
ful but one of the most useful features of the 
book is probably the excellent and thorough 
54-page bibliography extending over the 
1920-1960 period, with the full titles of the 
papers given. This bibliography is invaluable 
to all physicists engaged in this field of re- 
search. The main reservation in the opinion 
of this reader is that group theory was treated 
too briefly, maybe half-heartedly and that 
the elegance of group-theoretical methods has 
not been fully exploited. 

In summary, Slater’s two volumes on 
atomic structure are almost sure to become 
classical. While the first volume is quite 
suitable as a textbook, the second will be a 
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basic reference book. It is quite likely that 

many readers of Slater's books are anxious 

(like the present writer) to see the promised 

volumes III, IV, etc... . dealing with 
molecules and solids. 

HENRI AMAR 

Temple University 


PRINCIPLES OF METEORITICS, by E. L. Krinov. 
535 pages, illustrations, 54 X 84 in. Ox- 
ford, Pergamon Press, Ltd., 1960. Price, 
70s net. 


The advent of the space age has given new 
meaning to meteoritics. Today, by means 
of satellites and rocket probes, the scientist is 
consistently sampling space to uncover its 
content of particulate matter. 

In the past we were passive observers who 
could only study comets when they came into 
the sky, who could only see meteors when 
they flashed across our sky and who could 
examine meteroites only when they fell to 
Earth. As a result of our passive role the 
volume of information coming to us was 
strictly limited; yet even with these limita- 
tions there have been many in the field who 
have contributed significantly to our knowl- 
edge. Names like Olivier, Wylie, Ninninger 
and others come to mind when we discuss 
things that fall from the sky. However, these 
are American astronomers and researchers 
who have devoted their lives to this study. 
Others in other lands have also contributed 
significantly to this science. As Harrison 
Brown in his preface writes: ‘It [this book ] 
is admittedly nearly as deficient in reference 
to English literature on the subject as our 
works are deficient on reference to Russian 
literature.’ 

One of the strongest teams has been the 
Russians who have diligently pursued this 
science for a long time. Names like Kulik, 
Levin, Fesenkov and others in the Soviet 
Union are well known for their contributions. 
The author, like the others, is an authority 
in this field. 

Principles of Meteoritics is a well written, 
impressive treatise on meteoritics. The author 
has covered the various facets of this science 
in a lucid and brilliant fashion. 

The book begins in the traditional manner 
with the historical background of meteorites 
and then covers the motions and fall of these 
objects from the sky. In this section many 
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of the questions posed to astronomers such 
as the speed of meteors, sounds that ac- 
company them, brightness, sizes, etc., are 
answered. The mechanism involved in the 
burning out of a meteor is also covered though 
better explanations have appeared in English 
texts. Interesting tables going back to 1798 
are included indicating the fall of meteorites 
which have struck buildings or people. The 
most prominent of the meteor craters on the 
Earth are also described and discussed. In 
the case of the Tunguska meteorite the author 
indicates that the latest theory is that this 
meteorite was really the nucleus or head of a 
comet which struck the Earth to account for 
the odd effects felt on impact. 

One of the most significant sections of this 
book is on the morphology of meteorites. 
These studies of the dimension, form, surface 
relief, structure of the fusion crust and of the 
inner structure is most important in our under- 
standing of these objects. The author has 
performed a magnificent bit of writing in this 
section to account for the variations we find in 
meteorites. 

For the first time this reviewer has found 
an explanation as to why there are the vari- 
ations in the widmanstatten figures in meteor- 
ites. The author also gives a quick method 
of checking the authenticity of a meteorite 
by its nickel content. He says that while 
the nickel content of terrestrial iron does not 
exceed 2-3 per cent, no meteorite has been 
found with a nickel content of less than 4%. 
Some go as high as 30 per cent or even more. 

Following two sections on the chemistry of 
meteorites and the mineralogy of meteorites 
there is a section on the physical properties 
of meteorites. Here the author explores the 
specific gravity, melting points, magnetic 
properties and optical properties of meteorites. 

In two appendices are listed the meteorites 
which are in the USSR and a catalogue of 
rare meteorites. 

There is one omission of major significance. 
The author has not included in his book re- 
sults, current findings, on particulate matter 
in space which falls to the Earth. These 
results are rather definitive as they come from 
rocket probes and from satellites. Instead of 
the 6.6 tons falling to the Earth in a 24-hr. 
period, the latest results—which were avail- 
able to the author—indicate that between 
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1000 and 10,000 tons fall to the Earth in this 
period. 

Principles of Meteoritics is a brilliant work 
and promises to be a ranking reference book 
on meteoritics for a long time. 

I. M. Levitt 
The Fels Planetarium 


MECHANICS OF SOLIDS AND FLurps, by Robert 
R. Long. 156 pages, diagrams, 6 X 9 in. 
Englewood Cliffs, Prentice-Hall, Inc., 1961. 
Price: $9.00 (trade) ; $6.75 (classroom). 


Written for the instruction of undergradu- 
ate engineering students, this introduction to 
the analytical mechanics of continuous media 
makes consistent and effective use of vector 
and tensor methods, to which the introductory 
chapter is devoted. Their application to the 
solution of problems in stress, deformation 
and flow is illustrated in the two succeeding 
chapters. In accordance with the author’s 
intent to write an engineering book rather 
than one on mathematical physics, the fourth 
chapter discusses the mechanical properties 
of materials and points out the extent to which 
most of them differ from the two convenient 
abstractions, isotropic elastic solids and New- 
tonian fluids. 

The tightly yet clearly written fifth chapter 
derives the equations of elastic solids with 
the use of tensor notation in only thirteen 
pages, and provides an exhibition of its power. 
Chapter 6 illustrates the application of these 
equations to a selected group of problems in 
linear elasticity. Chapters 7 and 8 deal, re- 
spectively, with the equations of fluid me- 
chanics and their application to specific prob- 
lems in both statics and dynamics. 

The author has admirably succeeded in 
writing an introductory text on the mechanics 
of continua, which is rigorous and complete 
enough to serve as a foundation for the study 
of more extensive and specialized works. Yet 
the style is both lucid and informal, and clearly 
displays the author's skill in teaching a funda- 
mental engineering subject. There are many 
thought-provoking examples, and the refer- 
ences are not merely cited, but evaluated in 
their specific context, with the material dis- 
cussed in each chapter. 

CAPPELL 
The Franklin Institute Laboratories 
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LINEAR GRAPHS AND ELECTRICAL NETWORKS, 
by Sundaram Seshu and Myril B. Reed. 
315 pages, diagrams, 6 X 9 in. Reading 
(Mass.), Addison-Wesley Publishing Co., 
Inc., 1961. Price, $9.75. 


Introducing the theory of linear graphs to 
electrical engineers and demonstrating the 
power of these methods in solving network 
problems, are the aims of this text. 

The first five chapters are devoted to the 
basic theory of graphs in relation to their 
engineering applications. Considerable space 
is devoted to the matrices of a graph, and the 
relevant concepts are also treated in great 
detail. The last five chapters, which com- 
prise the majority of the book, treat appli- 
cations of graph theory. These deal with 
electric network theory, switching theory, 
communication networks, signal flow graphs, 
and logic. 

The book establishes a solid mathematical 
foundation for electric network theory, par- 
ticularly those formulations which are appli- 
cable to digital computers, and is based upon 
the results of several years’ research and 
teaching by the co-authors. 


CoBaLT: CHEMISTRY, METALLURGY AND USEs, 
edited by Roland S. Young. 424 pages, 
diagrams, 6 X 9 in. New York, Reinhold 
Publishing Corp., 1961. Price, $15.00. 


The chemist who has occasion to work with 
cobalt and requires accurate, up-to-date and 
complete information, and the chemist or 
metallurgist who wishes to browse through 
this field in the hope of adding to his general 
knowledge or acquiring ideas which may apply 
to his own sphere of work, are the ones to 
whom this book is directed. 

The latest information on the production 
and consumption of cobalt in all its forms, 
and a complete survey of the many new in- 
dustrial applications of cobalt that have been 
developed in the past decade, are found in 
this text. The properties and compounds of 
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this increasingly important metal are discussed 
by leading authorities in the field, in addition 
to the occurrence and extractive metallurgy 
of cobalt. 


MAGNETIC TAPE INSTRUMENTATION, by Gomer 
L. Davies. 263 pages, diagrams, 6 X 9 in. 
New York, McGraw-Hill Book Co., Inc., 
1961. Price, $8.50. 


This guide to the uses of magnetic tape 
reproduction includes information on tech- 
niques and equipment for reproduction of 
data (not voice and music), such as electrical 
signals or physical phenomena which can be 
represented by electrical signals. It is par- 
ticularly useful as a new technique for scien- 
tists and engineers, for it requires no proces- 
sing between recording and reproduction. 
The book is aimed at users of magnetic tape 
equipment, to give them a broader under- 
standing of this new tool. It is not a de- 
signer’s manual, but a user’s one, since appli- 
cations are stressed. 


SEPARATION PROCESSES IN PRACTICE, edited 
by Robert F. Chapman. 209 pages, dia- 
grams, 5 X 7} in. New York, Reinhold 
Publishing Corp., 1961. Price, $4.95. 


The latest information on both theory and 
practice of the various separation processes, 
based on fundamental processes, is presented 
in this small text. The newest theories are 
examined and presented in the first two chap- 
ters along with a review of the fundamentals. 
Typical examples of both the oldest and the 
newest techniques employed in the separation 
process are discussed in the remaining chap- 
ters. Described also is some of the latest 
equipment now in use in the field. This work 
is based on a symposium, “Experiment in 
Industry"’ which was jointly sponsored by 
the Philadelphia-Wilmington Section of the 
American Institute of Chemical Engineers 
and the University of Pennsylvania. The 
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engineer who desires to brush up on funda- 
mentals, broaden his practical knowledge of 
the field, and learn of new developments will 
find this book of valuable assistance. 


MoperN Fiicut Dynamics, by W. Richard 
Kolk. 288 pages, diagrams, 6 X 9 in. 
Englewood Cliffs (N. J.), Prentice-Hall, 
Inc., 1961. Price: $10.00 (trade); $7.50 
(text). 

This volume represents the first in the forth- 
coming Space Technology Series. Its treat- 
ment is as a text and reference in the analysis 
of flight-vehicle design for aeronautical and 
astronautical engineers or for anyone engaged 
in the design or fabrication of airplanes and 
missiles. Chapter One begins with Newton’s 
Second Law and Lagrange’s Equation as 
alternative ways of formulating the equations 
of motion. The author then introduces in 
the second chapter direction cosines and 
Eulerian angles for the purpose of orientation 
reference. These two introductory sections 
are then used as tools with which the dynamics 
of flight vehicles are analyzed. The author 
also, in dealing with aerodynamic forces and 
moments, has not reduced equations of motion 
to nondimensional status for purposes of 
generality as do most books on aeronautics, 
but has treated all motion as érue and seconds 
as constant quantities—not varying with 
speed, altitude, and wing loading. For pur- 
poses of clarity and east of continuity, the 
text retains the “‘body-fixed’’ axis system 
throughout, rather than using the ‘‘velocity- 
fixed’’ system in some problems. Important 
features include a complete analysis of the 
roll-coupling phenomenon, the detailed de- 
scription of manual and powered control 
systems for aircraft in Chapters 9 and 10, 
and in Chapter 12, the approximate solution 
for the gyroscopic motion of a spin-stabilized 
rocket. Many illustrative examples are used 
to emphasize the analyses, and in addition 
numerous line drawings, graphs, and tables 
have been included. 


TRANSISTOR Circuit ANALYsts, by Maurice 
V. Joyce and Kenneth V. Clarke. 461 
pages, diagrams, 6 X 9in. Reading (Mass.), 
Addison-Wesley Publishing Co., Inc., 1961. 
Price, $10.75. 


The basic methods of analysis involved in 
the understanding and design of junction 
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transistor circuitry are presented in this text. 
The transistor models used are restricted to 
a small number that are all easily interrelated. 
Initially, the simple low-frequency models 
are presented and are related to the under- 
lying physical processes that occur in the 
device. Then after the reader has completely 
familiarized himself with these simple models 
the complications, such as high-frequency 
effects, noise, or saturation effects, are intro- 
duced. Special features of this text include a 
unified approach to transistor biasing prob- 
lems, use of the current peak of the tunnel 
diode as a reference device, and an outline 
of the use of a minority charge density model 
in transistor switching and saturation calcu- 
lations. New material is also included on 
emitter-follower transient responser, video 
amplifier with maximum flatness and maxi- 
mum gain bandwidth products, and distrib- 
uted amplifiers. The book can best be used 
as a reference for the circuit designer and as 
a senior-graduate level text. 


ADVANCED DESIGN IN STRUCTURAL STEEL, by 
John E. Lothers. 583 pages, diagrams, 
6 X9in. Englewood Cliffs (N. J.), Pren- 
tice-Hall, Inc., 1961. Price: $15.00 (trade) ; 
$11.25 (text). 


Practical and relatively simple solutions to 
the numerous complex and critical problems 
faced by the modern steel design engineer are 
provided in this volume. Up-to-the-minute 
data on continuous and rigid frame steel build- 
ings and statically indeterminate trusses have 
been included. Wind and earthquake stresses, 
floor and stair design, and steel frame con- 
nections for multi-story buildings are also 
appraised. The author has included such 
topics as light gage structural steel sections, 
torsional buckling, allowable web sheer and 
tension members. Important information is 
also given regarding Lamella roofs and one- 
story rigid frame bents. Numerous illustra- 
tions to clarify factual material, and valuable 
bibliographies with sources for methods of 
advanced analysis have also been included. 
Meeting all the requirements and demands of 
specialists in structural engineering, this book 
will assist industrial concerns, technical in- 
stitutes, and government agencies, 


| 
2 4 | 
ck 
re 
| 
| 


434 


Puysics or SEMICONDUCTORS, by A. F. Ioffe. 
436 pages, diagrams, 6 X Yin. New York, 
Academic Press Inc., 1961. Price, $12.50. 


Translated from the Russian work of 1957, 
the English version of Physics of Semicon- 
ductors has been revised and supplemented 
by the author, who is Director of the Institute 
of Semiconductors for the U.S.S.R. Seven 
chapters cover thoroughly solid electrolytes, 
metals, electronic semiconductors, quantum 
theory of semiconductors, properties of the 
crystal lattice, methods of measuring the 
basic properties of semiconductors, and semi- 
conducting materials, A final chapter gives 
the author’s conclusions briefly and clearly. 


BooLEAN ALGEBRA AND ITs APPLICATIONS, by 
J. Eldon Whitesitt. 182 pages, diagrams, 
6 X9in. Reading (Mass.), Addison-Wes- 
ley Publishing Co., Inc., 1961. Price, 
$6.75. 


Designed as an introduction to Boolean 
algebra for readers with a limited formal 
mathematical background, this text can be 
used in a one-semester course. It is also 
suitable as background material for specialized 
courses in the three major applications of 
Boolean algebra—symbolic logic, probability 
theory, and circuit design. The seven chap- 
ters cover algebra of sets, the disjunctive and 
conjunctive forms of Boolean algebra, sym- 
bolic logic and the algebra of propositions, 
switching algebra, relay circuits and control 
problems, circuits for arithmetic computation, 
and probability in finite sample spaces. 


INTERFEROMETRY. National Physical Labor- 
atory Symposium No. 11. 471 pages, il- 
lustrations, 64 X 94in. London, Her Maj- 
esty’s Stationery Office, 1960. Price, $5.60 
postpaid. (U. S. agent, British Informa- 
tion Services. ) 


Interferometry is the proceedings of a sym- 
posium held at the National Physical Labora- 


tory in June, 1959. Nineteen papers and 
the resulting discussions comprise the text. 
Since all the papers were circulated in ad- 
vance, the discussions form a valuable part 
of the proceedings. Primarily for physicists, 
the symposium covered recent developments 
in interferometry, particularly as related to 
metrology, optics, radio and radiation. 
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PERCHLORATES: THEIR PROPERTIES, MANU- 
FACTURE AND UsEs, edited by Joseph C. 
Schumacher. ACS Monograph No. 146. 
257 pages, diagrams, 6 X 9in. New York, 
Reinhold Publishing Corp., 1960. Price, 
$8.75. 


All the current knowledge of the technology 
of perchlorates and perchloric acid has been 
brought together in this new ACS monograph, 
by leading authorities in the field. They have 
traced the commercial development of the 
perchlorate industry on a world-wide scale 
and reviewed the properties and manufacture 
of the compounds, discussing their safe 
handling and examining their various uses. 
An entire chapter is devoted to the application 
of perchlorates in the explosive and propellant 
industries. This material is particularly timely 
due to the growing interest in the compounds 
as oxidizers in propellants for the military 
space program. This book, aside from its 
immediate value as an accurate reference 
for present day problems encountered in the 
manufacture and use of perchlorates and 
perchloric acid, provides the background for 
future research and development in this ex- 
panding field. 


Linear Circuits, by R. E. Scott. Part I, 
TrmE-DomatIn ANALYysts, 510 pages, dia- 
grams, 6 X9 in. Price, $6.75. Part II, 
FREQUENCY-DoMAIN ANALYsiIs, 418 pages, 
diagrams, 6 X 9 in. Price, $6.75. Read- 
ing (Mass.), Addison-Wesley Publishing 
Co., Inc., 1960. 


This introduction to linear circuits for 
students of electrical engineering presents 
the methods of linear circuits as a solid founda- 
tion upon which the complete structure of 
electrical engineering can be built. The 
author assumes in this work that the student 
has had an introductory course in physics 
and is studying calculus concurrently. Part 
I, containing fourteen chapters, treats the 
basic theory, and topics such as singularity 
functions, response functions in the time 
domain (that is, step and impulse responses), 
differential equations, and convolution theory. 
Part II, with nine chapters, treats the network 
system functions in the frequency-domain 
(that is, responses to steady-state sinusoids), 
and Laplace transforms and power density 
spectra for more general input functions. 
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The text also covers topics seldom found in 
books at this level, such as signal-flow graphs, 
relaxation methods, dummy variables, re- 
sponse functions in the time and frequency 
domain, and s-plane plots. 


INTRODUCTION TO THE THEORY OF IONIZED 
Gases, by J. L. Delcroix. Translated by 
Melville Clark, Jr. and David J. BenDaniel. 
149 pages, diagrams, 5} X 8in. New York, 
Interscience Publishers, Inc., 1960. Price: 
$4.50 (cloth) ; $2.50 (paper). 

Translated from the 1959 French edition, 
this text approaches the theory of ionized 
gases through a microscopic description of the 
state of the gas. This complements Spitzer's 
well known work, using macroscopic equa- 
tions. The material has been tested in 
courses given by the author at the Faculté 
des Sciences (Paris) and the Institut des 
Sciences et Techniques Nucléaires (Saclay). 
Twelve chapters deal with such topics as 
elastic collisions, the kinetic theory of gases, 
the theory of free electron gases, weakly 
ionized gases, and plasmas. Four short ap- 
pendices give some of the mathematical deri- 
vations. This is one of the few available 
works in the new field of plasma dynamics. 


52 FourrER TRANSFORMS, by R. R. Goldberg. 
76 pages, 54 X 8$ in. New York, Cam- 
bridge University Press, 1961. Price, $3.75. 


This study provides a background in clas- 
sical Fourier transform theory, in a form 
which lends itself to an understanding of ab- 
stract harmonic analysis. The introductory 
chapter contains all the theorems on integra- 
tion which are used in the remaining four 
chapters dealing with: the Fourier transform 
on L'; the Fourier transform on L?; generali- 
zations of Wiener’s theorem; and Bochner’s 
theorem. An appendix points out the rela- 
tionship to an arbitrary locally compact 
abelian group. The reader is assumed to 
have a knowledge of Lebesgue and Riemann- 
Stieltjes integration. 


INDUSTRIAL ORGANIC NITROGEN COMPOUNDS, 
by Melvin J. Astle. ACS Monograph No. 
150. 392 pages, diagrams, 6 X 9in. New 
York, Reinhold Publishing Corp., 1961. 
Price, $14.00. 


This monograph summarizes the chemistry 
of organic nitrogen compounds, with chapters 
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on: the aliphatic amines; arylamines; hetero- 
cyclic amines; hydrazine, azo compounds, 
diazonium salts and oximes; nitrides, amides 
and aminoacids; isocyanates, ureas, thioureas; 
aromatic nitro and nitroso compounds; anda- 
liphatic nitro compounds. Patent literature 
receives a great deal of attention, so that 
industrial as well as theoretical developments 
are included. The simple heterocyclic nitro- 
gen ring is discussed in detail, as background 
for an understanding of the nitrogen com- 
pounds. Some topics, such as amino acids 
and proteins, alkaloids, and heterocyclic 
nitrogen compounds, do not receive as much 
attention as the others, since they have been 
covered thoroughly in other works. 


HyYPERSTATIC STRUCTURES. VOL. 2: WORKED 
EXAMPLES AND EXAMPLES FOR SOLUTION, 
by J. A. L. Matheson and A. J. Francis. 
282 pages, diagrams 6 X 94 in. New 
York, Academic Press Inc., 1960. Price, 
$11.00. 


This work has been published in two com- 
plementary volumes, the second volume con- 
taining worked examples and problems apply- 
ing the principles and methods discussed in 
Vol. 1. The subject of hyperstatic structures 
(a structure is said to be “statically indeter- 
minate” or ‘‘hyperstatic’’ when it possesses 
more members or is supported by more re- 
active restraints than are strictly necessary 
for stability) is treated as an organic and 
coherent whole by the authors. The relation- 
ship between the various energy theorems is 
demonstrated by a consideration of non-linear 
structures; then the theorems and methods 
applying only to linear structures are treated 
in their logical order, with special attention 
being paid to the broad distinction between 
the “equilibrium’’ and ‘‘compatibility’’ ap- 
proaches. As an introduction to the analysis 
of structures on the digital computer, a 
chapter on matrix methods is included. 


ELEMENTS OF STATISTICAL ]NFERENCE, by 


David V. Huntsberger. 291 pages, dia- 
grams, 64 X 9} in. Boston, Allyn and 
Bacon, Inc., 1961. Price, $7.25. 


This book could serve as a text for a general 
course for those who want to know “what 
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statistics is all about.”” The techniques and 
methods which are presented are those which 
are generally applicable, and methods which 
are of interest in but a few areas of application 
have not been included. Some of the ex- 
amples and problems set forth in the book are 
artificial, some have been abridged from real 
sets of data, and others are the unabridged 
results of samples and experiments. 

A knowledge of arithmetic and some ex- 
posure to algebra are sufficient for the reader, 
since the mathematical prerequisites have 
been kept to a minimum. Those who want 
to be familiar with at least the language and 
fundamental concepts of statistics in order to 
achieve additional competence in their chosen 
fields would find this work beneficial. In 
summary, the author intended this book as 
an introduction to some of the basic concepts 
and techniques of statistical inference that 
play an important role in the solution of 
problems that arise in many spheres of human 
activity, and he hopes to give the reader 
some understanding of the statistical approach 
to problems faced by the administrator, the 
research worker, and other members of our 
increasingly complex society. 


Book Notes 


(J. F. 1 


NUCLEAR ENGINEERING, by Gilbert Cahen 
and Pierre Treille, translated by Gilbert 
B. Melese. 394 pages, illustrations, 6} X 8} 
in. Boston, Allyn and Bacon, Inc., 1961. 
Price, $8.50. 


Published originally as the highly successful 
“Précis d’Energie Nucléaire’ by the two 
authors listed above, two of France’s leading 
nuclear engineers, this book has been trans- 
lated and brought up-to-date by Dr. Melese 
to include the latest available information 
plus pertinent data from the 1958 Geneva 
International Conference. 

Those elements of modern physics which 
are significant in nuclear engineering are 
discussed in the first two chapters, followed 
in succeeding chapters by discussion on 
nuclear reactor theory, nuclear reactor ma- 
terials, and French, British and Russian re- 
actors. The book offers conversion tables, 
nuclear data tables, charts and graphs, and 
concludes with an examination of the hazards 
as well as the beneficial uses of radiation and 
includes a discussion of radiation shielding. 
The original French version of this work has 
been reviewed as “The best introduction to 
the immense domain of nuclear energy,’’ and 
has been widely used in Europe with great 
success. 


PUBLICATIONS RECEIVED 


SCIENTIFIC RussIAN, by James W. Perry. Second edition, 565 pages, 6 X 9 in. 
Price, $9.50. 


Interscience Publishers, Inc., 1961. 


New York, 


CHARACTERISATION OF ORGANIC CoMPoUNDS, by F. Wild. Second edition, 306 pages, 53 X 84 
in. New York, Cambridge University Press, 1961. Price, $2.95 (paper). 


Second edition, 520 pages, 6 X 9 in. 
Price: $12.00 (trade) ; $9.00 (text). 


ADVANCED CaLcu.us, by David V. Widder. Englewood 


Cliffs (N. J.), Prentice-Hall, Inc., 1961. 


THERMODYNAMICS, by G. N. Lewis and M. Randall, revised by Kenneth S. Pitzer and Leo 
Brewer. 723 pages, diagrams, 6 X 9in. New York, McGraw-Hill Book Co., Inc., 1961. 
Price, $12.50. 


DICTIONARY OF MECHANICAL ENGINEERING, by Alfred Del Vecchio. 316 pages, 6 X 9} in. 
New York, Philosophical Library, 1961. Price, $6.00. 


DatLy RIVER StaGEs, VoL. LIV, 1958. 174 pages, 8} X 10}in. Washington, U.S. Dept. of 
Commerce, 1960. For sale by the Superintendent of Documents, U. S. Gov’t Printing 
Office, $1.25 (paper). 
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NATIONAL BUREAU OF STANDARDS NEWS 


X-RAY BARRIER DETERMINATIONS 


Measurements made by the National Bureau of Standards, in a 
program sponsored by the Atomic Energy Commission, indicate that 
present methods of determining barrier requirements for X-ray installa- 
tions may result in walls having protective barriers that are several 
half-value layers thicker than actually required.! This study, con- 
ducted by S. W. Smith and J. R. Brooks, with the cooperation of over 
20 hospitals and ‘clinics, may lead to savings in the cost of many future 
installations. 

The usual procedure for the calculation of required barrier thickness 
is to assume near-maximum values for the various factors affecting the 
exposure of the barriers if more precise knowledge of them is not known. 
These factors include such variables as workload, fraction of time the 
radiation is directed toward the barrier, effect of coning, attenuation of 
the patient, and others. In actual practice, one or more of these factors 
may be considerably less than the maximum values assumed for busy 
installations and the barrier thickness required for personnel protection 
may thus be less than that indicated by the calculations based on these 
assumptions. 

In this study, the amount of radiation reaching various regions of 
the X-ray rooms was determined through the use of pocket-type ioniza- 
tion chambers. As many as 70 of these instruments were attached to 
selected positions on the ceiling and walls of each room. Readings were 
taken once or twice a week, depending on the exposure, for a period of 
from one to two months. The penetrating quality of the radiation was 
determined from attenuation curves obtained from readings of instru- 
ments that were surrounded with varying thicknesses of aluminum 
or lead. 

The work load in milliampere-minutes of each installation was auto- 
matically determined by a charge-integrating device connected to the 
X-ray control. From the data obtained it seems reasonable to assume, 
when no specific information is available, a workload of 1000 ma-min/ 
week for a busy diagnostic X-ray room. 

The amounts of lead shielding required to reduce exposure to the 
desired level were calculated from the collected data for those rooms 
having maximum exposure and use. These thicknesses were generally 


1 For further technical details see, ‘Measurement of Radiation Exposure at the Walls of 
Medical X-Ray Rooms for Determining Protective Barrier Requirements,”’ by S. W. Smith 
and J. R. Brooks, Atompraxis, March, 1960. 
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considerably lower than values arrived at using the usual procedure, 
indicating substantial safety factors. 

The results of this study, combined with data from other experi- 
ments, may eventually result in a reduction in the amount of shielding 
specified for X-ray rooms by use of more realistic values for the various 
factors used in computing the barriers. 


YOU CAN ADVANCE SCIENCE EDUCATION 


Today, more than ever before in its 137-year history, there is vital need for 
The Franklin Institute to effectively promote education in science and technology. 
It is imperative that we meet this challenge by providing adequate educational op- 
portunities in these fields. This requires vision, objective planning, and money. 
We have more than enough of the first two requisites, but far too little of the third. 

Our programs are aimed at professional scientists and industry, as well as the lay 
public and young people seeking inspiration and guidance in choosing a career. 
The Institute’s educational programs are impressive. They begin with students in 
the early grades of our elementary schools and continue throughout an individual's 
professional or industrial life. With more funds at our disposal, the scope and 
vigor of these activities could be greatly increased and increasingly effective. 

The Franklin Institute is not richly endowed. It is a non-profit organization, 
depending for encouragement and support on an understanding public. Capable 
and conservative management assures wise administration of all funds. 

Your gift or bequest, large or small, will be deeply appreciated and will be used 
effectively to broaden the Institute’s educational usefulness. There is a warm 
satisfaction in giving financial support to an organization that has pioneered in, and 
is dedicated to, the advancement of science and technology. 

When property is transferred, title should be in the name of The Franklin 
Institute of the State of Pennsylvania for the Promotion of the Mechanic Arts. 

The Secretary of The Franklin Institute will gladly furnish you with additional 
information. Write to him at The Franklin Institute, Benjamin Franklin Parkway 
at Twentieth Street, Philadelphia 3, Pennsylvania. 
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CURRENT TOPICS 


“Phantom Holes’ in Computer 
Tape.—A problem caused by oil spots 
which, in effect, became “phantom 
holes’”’ in paper tape for instructing 
computers, and led to faulty data- 
processing results, has been solved by 
the Radio Corporation of America. 

Solution of this problem—common 
throughout the computer industry— 
has been achieved by RCA through 
development of a high-speed tape 
reader that knows a genuine perfora- 
tion when it sees one. 

Much of the paper tape employed 
in present-day electronic data proces- 
sing is oil-impregnated and small areas 
of oil concentration sometimes occur. 
Oil spots also may be caused by 
lubricants in connection with high- 
speed operation of data-processing 
equipment. 

Photoelectric scanning principles 
and improved circuitry have been in- 


corporated in the design of RCA’s . 


Model 322 Tape Reader, enabling 
it to ignore everything but true 
perforations. 

In testing the selectivity of the new 
tape reader, RCA engineers deliber- 
ately sprinkled perforated paper tapes 
with oil and then ran them through 
the tape reader. The results were 
uniformly successful, he said. 

There are two models of the 
new RCA paper tape equipment—one 
serves as a tape reader only, while 
the other can read as well as punch 
the data-signifying holes. The reader 
can scan perforated tape at rates up 
to 1,000 data characters per second, 
handling 5, 6 or 7-level formats. Ac- 
curacy in the reading operation is as- 
sured by a parity check. 

Resembling an undersized secretary 


desk, the Model 322 Paper Tape 
Reader/Puncher stands 60 in. high 
and measures only 34 in. by 24 in. 
at the base. 


New Atomic Clock.—An atomic 
clock, so accurate that its maximum 
error would not exceed one second in 
1271 years, has been developed for 
the Air Force. The airborne device 
which uses atomic energy as a power 
source is technically referred to as an 
Airborne Atomic Frequency Standard. 
It should become operational for use 
in Air Force missiles and aircraft next 
year. 

The 623-lb. airborne atomic fre- 
quency standard will replace the nu- 
merous crystal oscillators, used as 
frequency or time standards, now re- 
quired to calibrate airborne communi- 
cations, navigation, guidance, fire con- 
trol, computers and timing devices. 

First of its kind to be developed, 
the unit has been simplified to the 
extent that only a simple on-off switch 
is required to operate it. Present 
oscillators require frequency calibra- 
tion, often just before takeoff or during 
flight. The new atomic standard does 
not require calibration, acting as its 
own primary standard. Operator of 
the unit needs no special schooling 
such as is required on present labora- 
tory units. 

The small clock not only will be 
airborne, but can also replace the 
present-day atomic laboratory units 
weighing over 600 pounds. No radia- 
tion hazard is involved in its operation. 

The airborne atomic frequency 
standard is accurate to one part in 
10 billion and is stable to one part in 
25 billion. 
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Combining the latest techniques in 
solid state devices, the atomic fre- 
quency standard was developed to air- 
borne specifications by the National 
Company, Inc., Malden, Mass., under 
the guidance of the Communication 
Laboratory, Wright Air Development 
Division, Air Research and Develop- 
ment Command, Wright-Patterson 
AFB, Dayton, Ohio. 


Telemetry System for Space.— 
General Electric’s Missile and Space 
Vehicle Department in Philadelphia 
has completed an experimental version 
of a new telemetry system that is 
expected to be able to beam signals 
through space about three times as 
far as any telemetry system yet flown, 
including Pioneer V, United States 
solar satellite which last year exceeded 
all records by sending signals to Earth 
from millions of miles out. The tele- 
metry system is one of a family of 
communication systems called Syn- 
chrolink that will transmit data the 
same distance as systems now in use 
with about one-tenth the power re- 
quirements. The estimated ranges are 
based on the inherent efficiencies of the 
system, which is described as ‘“‘the 
most efficient practically realizable 
with today’s technical know-how.” 

Synchrolink is a digital system, 
which gives it a high degree of data 
accuracy. It is classified as Pulse 
Code Modulation with Phase Shift 
Keying (PCM/PS). An optional fea- 
ture is an error correction code which 
can be used to correct mistakes caused 
by noise interference. 

Other reasons Synchrolink is at- 
tractive are: 

1. Relatively minor modifications 
would be needed to convert existing 
ground stations to receive Synchrolink 
signals. 

2. Transistors are used in their most 
reliable mode: as switches. 
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3. It is flexible: channels and com- 
binations of channels could be altered 
to accommodate measurement changes, 
even in flight. 

4. Repeated conversion, or trans- 
mission of data, or long storage, causes 
no loss in accuracy. 

Ina typical application, information 
such as the space vehicle performance 
and temperature environment is gath- 
ered by sensors located throughout 
the vehicle. Samples of information 
direct from the sensors, or via pre- 
transmission analyzers, are coded into 
a series of positive or negative pulses. 
The Synchrolink system then beams 
a radio frequency signal that is 
switched 180 degrees in accordance 
with a simple coded pulse series repre- 
senting the space flight data. The 
resulting double side band-suppressed 
carrier radio frequency signal is re- 
ceived by a synchronous receiver fol- 
lowed by a correlation detector. Use 
of the threshold-free synchronous re- 
ceiver permits detection of the signals 
at extremely low signal-to-noise ratios, 
about —10 db. 

The detected output is recorded on 
tape and run through decoding and 
processing equipment most of which 
is standard. 

Components contributing to the 
unique properties of the Synchrolink 
system are: Error Correction Encoder, 
Modulator for Phase Shift Keying, 
Synchronous Receiver, Correlation 
Detector, and Error Correction 
Decoder. 


Classroom Alpha Particle Counter. 
—A new Open Air Alpha Counter for 
demonstrations in the high school and 
junior college classroom has been 
announced by Central Scientific, 
Chicago. 

The apparatus is able to detect 
alpha particles as they pass through 
the instrument causing a dielectric 
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breakdown of air under a high po- 
tential gradient. This breakdown, as 
indicated by a small spark, occurs 
whenever a particle ionizes a path in 
the air between the Counter’s nickel- 
chromium grid and detector plate. 
These two are maintained at about 
5000-volt d-c. potential difference by 
an accessory power supply. The 
source of alpha radiation (radium 
salts) is mounted on an adjustable 
nylon rod that can be moved back and 
forth in relation to the counter grid. 
By varying the d.-c. potential or the 
position of the source, the rate of 
counts can be closely controlled. When 
operated with a scaler, an accurate 
count of particles can be made over a 
pre-determined time interval. 

A clear plastic cover protects the 
radioactive source and counting grid 
from student contact and dust during 
storage. The Counter is mounted on 
a vertical rod in a stable Harrington 


base. The instrument requires little 
maintenance other than periodic clean- 
ing of the grid wires with carbon 
tetrachloride. 

The counter operates at a maximum 
potential of 5000 volts d.c. at a maxi- 


mum of 1 milliampere. A variable 
5000-volt d-c. source is also available 
from Cenco. 


Moisture Detector for Missile Sys- 
tems.—A moisture detector has been 
developed by the U. S. Army Engineer 
Research and Development Labora- 
tories, Fort Belvoir, Va., for use at 
missile installations and other areas 
where undetected moisture might 
cause malfunction of vital instruments 
and equipment. 

Called an “electrolytic hydrom- 
eter,’ the instrument is capable of 
detecting moisture in high pressure 
air up to 6000 psi. and also of measur- 
ing moisture content down to a free 
air dew point of — 100° F. It weighs 
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35 lb., and is 11 in. high, 12 in. wide, 
and 15 in. long. 

The instrument consists basically 
of an electrolytic cell, an ammeter 
and a direct current power source 
with associated pneumatic and elec- 
trical controls for efficient operation. 

In use, the electrolytic type mois- 
ture analyzer absorbs and simultane- 
ously electrolyzes all moisture from 
a sample gas stream. The absorption 
and electrolysis are accomplished in a 
small diameter tube containing a pair 
of closely spaced platinum wires 
wound as a double helix on its inner 
surface and coated with a thin film of 
phosphorous pentoxide, a material 
with a strong affinity for moisture. 

The equipment was built by Con- 
solidated Electrodynamics Corp. and 
Beckman Instruments Inc., under 
contracts with the Laboratories. An- 
alyzers now in use are with the Red- 
stone missile systems. 


Electronic Tape Converters.— 
Minneapolis-Honeywell’s Electronic 
Data Processing division has an- 
nounced the development of electronic 
converters that automatically trans- 
late records on magnetic tape from 
IBM and Univac language into 
Honeywell language. The IBM con- 
verter has been completed and is in 
operation at the division’s service 
bureau at Wellesley, Mass., where 
the Univac converter will soon be in 
operation. 

For the present, Honeywell plans 
to use the new machines on a service 
basis to assist those users of Honeywell 
EDP systems who wish to convert 
non-Honeywell tapes into Honeywell 
tapes. 

The translations proceed at what- 
ever speed the non-Honeywell tapes 
read the recorded data. The trans- 
lations often result in substantial 
tape savings, with the ratio running 
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as high as two and one-half non- 
Honeywell tapes to one Honeywell 
tape. 

The equipment operates as follows: 

Univac and IBM tape units read 
the information recorded on their 
tapes into the Honeywell converters 
(the Model 831 for Univac and the 
Model 832 for IBM). This equip- 
ment converts the non-Honeywell 
signals into Honeywell words and pas- 
ses them through a standard Honey- 
well tape control unit into the central 
processor of a Honeywell 800. The 
Honeywell 800, assisted by one of 
several specially-written programs, 
translates the IBM or Univac code 
into Honeywell code and records the 
output on Honeywell magnetic tapes. 


Magnetic Mathboard Simplifies In- 
struction.— Mathematics can be more 
effectively presented to students with 


the impact gained from visualization 
with the new Magnetic Mathboard 
now available from Atomic Labora- 
tories, Inc. 

The Mathboard is a 3 ft. by 4 ft. 
magnetic blackboard that can be set 
on its own easel or on the chalk tray 
of a regular classroom blackboard. 
By means of white and black plastic 
geometric figures magnetized to the 
board, scores of mathematical func- 
tions can be visually illustrated by 
moving them around the board. Mov- 
able angle pieces, for example, show 
that the sum of a triangle’s three 
angles equals 180 degrees. 

In arithmetic, the Magnetic Math- 
board illustrates slide rule operations 
of addition and subtraction, frac- 
tions and their manipulation; in 
algebra, it shows algebraic addition 
and subtraction, monomial multipli- 
cations, factoring the difference of two 
squares, graph work, arithmetic pro- 
gression; in geometry, it illustrates 
the sum of angles in geometric shapes 
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and other geometric functions; and 
in trigonometry, it shows multiplica- 
tion and division using logarithms. 
A 26-page, fully illustrated manual 
provides a complete selection of 
demonstrations. 

In addition to holding the various 
magnetized plastic parts, geometric 
shapes, numbers, letters, etc., the 
board functions as a regular black- 
board and can be written on with 
chalk. The board comes with a 
regulation blackboard finish and is 
framed in aluminum trim. 


Data-Processing System for Bank- 
ing.—Burroughs Corporation has an- 
nounced a new electronic system for 
speeding and controlling the inter- 
change of more than 50 million checks 
every day among U.S. banks. Called 
the Burroughs B270 Electronic Data 
Processing System, it represents 
another major advance in bank auto- 
mation equipment developed to curb 
rising operating costs and to alleviate 
a shortage of trained personnel. 

Designed to provide automatic 
proof and transit operations, the B270 
will also convert ‘“‘on-us” items to 
magnetic tapes for subsequent com- 
puter processing and may be used for 
automatic deposit analysis as well as 
account reconciliation services. 

The system's versatility will enable 
banks to expand their correspondent 
services and improve item control for 
either conventional bookkeeping or 
electronic data-processing systems, ac- 
cording to the company. 

With ever growing volumes of 
checks flowing through their proof and 
transit departments, many banks are 
finding it increasingly difficult to 
maintain adequate staffs of well quali- 
fied proof clerks. 

It takes several months to train a 
skilled operator of conventional proof 
machines, with much of this time 
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devoted to teaching complicated rout- 
ing patterns. Automatic routing on 
the B270 eliminates this problem. 

Furthermore, the system offers an 
answer to the profit squeeze caused by 
rising proof and transit costs through 
reduced item handling and greater 
over-all efficiency and accuracy. 

The B270 combines magnetic ink 
character recognition (MICR)—bank- 
ing’s ‘common machine language’”’— 
with high speed document sorting, 
multiple-list printing and solid-state 
electronics. This results in broad item 
distributions with a minimum of 
handling. 

The system utilizes the world’s 
fastest electronic sorter-reader, which 
reads checks and deposit tickets en- 
coded with MICR at speeds up to 
1560 items per minute, for input to 
the system. The sorter is also used 
to accomplish the physical sorting of 
all items. 


Other units include an internally- 
programmed central processor with 
magnetic core memory, two ultra-high 


speed multiple-tape listers and a 
punch card reader for program loading 
and entry of miscellaneous informa- 
tion into the system. 

Magnetic tape units are added to 
the system to permit conversion of 
“on-us’”’ information for subsequent 
use in demand deposit data processing. 
These units also facilitate the use of 
the B270 as an off-line converter for 
various data processing systems. 

Improved item control stems from 
the system’s ability to produce master 
listings and detailed cash letters in 
which each item is identified with 
either the payor bank or the individual 
payor. This extended audit trail is 
obtained automatically, with no re- 
duction in handling speed. 

Application concepts of the B270 
have been proven in Burroughs in- 
stallations of magnetic document pro- 
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cessing systems now in operaion at 
Michigan National Bank of Lansing, 
First Pennsylvania Banking & Trust 
Co. of Philadelphia and the Chicago 
Federal Reserve Bank. 


Epoxy Silver Solder.—A silver con- 
ductive epoxy with a resistivity ap- 
proaching the resistance of metals— 
between 0.01 ohms-centimeter and 
0.0001 ohms-centimeter—has been de- 
veloped by Epoxy Products, Inc., a 
division of Joseph Waldman & Sons, 
137 Coit Street, Irvington, N. J. 
Epoxy solder has two principal ad- 
vantages over conventional solder or 
brazing. It can be used at low tem- 
perature for components which are 
heat sensitive and it has no flux or 
residue to contaminate sensitive elec- 
tronic components. 

Epoxy solder is available in two 
paste forms—a one component heat 
curing (as low as 125° C.) paste and 
a two component room temperature 
cure paste. It is ideal for field repair 
work, particularly on crowded printed 
circuits where solder temperatures 
might damage heat-sensitive compo- 
nents. Its bond strength is superior 
to usual metal solders and it can be 
used to make a conductive bond with 
practically any material. Shear 
strength of a steel-to-steel bond is 
3200 psi. Aluminum has approxi- 
mately the same shear strength but 
is difficult to measure because the 
aluminum begins to tear before the 
bond fails. 

When the bond hardens it can be 
planted with conventional plating 
techniques. Soldering can then be 
performed on the plated surface if 
desired. 

Unlike conventional soldering and 
brazing, epoxy solder requires no flux 
and leaves no residue which might 
contaminate sensitive electronic com- 
ponents. Typical applications in- 
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clude bonding electronic tube caps to 
the glass tube, bonding the tantalum 
slug to the can of a tantalum capacitor, 
connections made to _heat-sensitive 
diodes and transistors and grounding 
aluminum capacitors to steel chassis. 

Conductive epoxies are intended 
primarily for replacing soldering and 
brazing operations which require high 
temperatures. When used with dis- 
similar metals, epoxy solder offers a 
bonding technique where no type of 
soldering or brazing would be possible. 


Solid-State Crystal Oscillator.—A 
miniature, all solid-state voltage-con- 
trolled crystal oscillator (VCXO), the 
first standard model of a new line of 
similar products, was shown by Itek 
Electro-Products Company at the 
1961 IRE Show, held in New York 
recently. 

The Itek development gives radio 
engineers, for the first time, a crystal- 
controlled FM oscillator which can be 
directly frequency-modulated over a 
range as great as + 0.2 per cent of 
center frequency, with excellent lin- 
earity and with virtually no degrada- 
tion of crystal frequency stability over 
a wide range of temperatures. The 
units are expected to find wide appli- 
cation for FM generation, automatic 
frequency control and other opera- 
tions in communications, telemetry, 
and guided missile systems. 

VXCO Model 10M-WA has a center 
frequency of 10.7 mc and a deviation 
of + 20 kc; linear to within + 200 
cps. Measured temperature drift is 
less than 1 kc. from — 40 to + 65° C. 
Output power is 5 mw.; size of the 
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transistorized unit is only ? in. X2 in. 
X 23 in. over-all. 

The standard Model 10M-WA is 
designed with a center frequency of 
10.7 mc. to be compatible with Itek 
Electro-Products’ 10.7 mc. Series of 
high-selectivity Hermes crystal filters 
and matching crystal discriminators. 
The 10M-WA is available in limited 
quantities for immediate delivery; 
other VCXO’s, in the center frequency 
range 10 kc. to 30 mc. can be built to 
special order to meet both commercial 
and military system requirements. 
Units can also be supplied housed in 
temperature-controlled ovens which 
minimize temperature drift and in- 
crease frequency stability. 

The VCXO circuit incorporates a 
unique embodiment of crystal reso- 
nator and voltage-variable reactance 
to permit direct frequency modula- 
tion of the oscillator. This technical 
achievement derives from a thorough 
mathematical analysis, by Itek Elec- 
tro-Products’ engineers, of the prop- 
erties of piezoelectric crystals used 
with variable reactance networks of 
various types. The device is patented. 
It is claimed to be no more complex 
in size and number of components 
than conventional amplitude modu- 
lating systems, permitting the ad- 
vantages of frequency modulation in 
many crystal-controlled applications 
where it had not been feasible before. 
This is especially true of portable 
equipment, where the conventional 
FM technique of multi-stage phase 
modulation synthesis requires con- 
siderably more complex circuitry and 
larger equipments. 
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STAFF POSITIONS 


... are available for qualified scientists and engi- 
neers in the Laboratories of The Franklin Insti- 
tute. These Laboratories comprise one of the 
country’s finest independent research organiza- 
tions, and provide R&D services to both in- 
dustry and government, in a variety of fields: 


chemistry, solid state sciences, applied physics, 
operations analysis, engineering psychology, 
and chemical, electrical, mechanical and nu- 
clear engineering. 


There are openings, in groups working on chal- 


lenging problems in both research and develop- 
ment, for men with proven abilities in any of 
chese fields. 


THE FRANKLIN INSTITUTE LABORATORIES 
FOR RESEARCH AND DEVELOPMENT 


20th Street and Benjamin Franklin Parkway 
Philadelphia 3, Pennsylvania . .. LOcust 4-3600 


Send complete resumé to 
Mr. John E. Christ, Director of Personnel 


Or write for a brochure describing activities 
of the Laboratories 


ip 
Petes: 
= 
‘ 
: 


BUSINESS 


MATH ALY 


ENGINEERING 


Now in the fourth year of operation, The Franklin Institute Com- 
puting Center has shared its technical know-how, modern electronic 
equipment and experienced personnel with hundreds of progressive 
industries and government agencies across the nation. 


We offer the services of creative people, highly ski!led in their respec- 
tive fields and ably trained in the application of these skills to the ever 
expanding area of electronic computers and data processing systems. 


This staff is now available for analysis, system design, programming 
or coding of projects of unlimited scope or context. Input to our large 
scale computer and peripheral equipment is acceptable in any form. 
Results are provided on cards, plastic or metallic tape, and in com- 
pletely edited printed copy. Our extensive library of automatic coding, 
engineering, data processing and mathematical routines is available to 
all users, and machine time is provided with or without the services 
of programming personnel. 


THE FRANKLIN INSTITUTE 
Computing Center 


20th Street and Benjamin Franklin Parkway 
Philadelphia 3, Pa. LOcust 4-3600 
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